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ASTP EXPERIMENT SUPPORT DATA PROCESSING 


By Richard K. Osburn, 

Edward L. Barnett, Harry L. Moore, 
John B. Moore and Jimmie R. Ball 


1.0 INTRODUCTION 


This internal note documents in detail the activities associated with the 
generation of ASTP experiment support data in the areas of- spacecraft ephemeris 
and orientation and instrument pointing and field-of-yiew. The majority of this 
work is being performed for the Mathematical Physics Branch (MPB) by TRW engineers 
under contract NAS 9-13834, Task JSC/TRW 308. It is intended that this document 
represent a cradle-to-grave chronicle of these activities. To satisfy this intent 
while facilitating the ready dissemination of information, the document is being 
published twice. The first publication, scheduled for release prior to ASTP lift- 
off, includes all preflight phases of the experiment support activity in addition 
to those appendixes that do not pertain to any mission-specific data. The second 
publication will provide any required updates to the original documentation and 
will add all mission-specific data, including documentation of all postflight 
data processing activities and data archiving information. 

It should be noted that the responsibilities of MPB and Task 308 are limited 
to the generation of spacecraft ephemeris, instrument pointing and field-of-view, 
and related experiment support data. The task is not involved in the retrieval 
of CSM orientation data, experiment sensor data, or any other such data from raw 
telemetry. In fact, the orientation data must be recovered by the IDSD of DSAD 
and delivered to MPB before generation of experiment support products can begin. 



2.0 EXPERIMENT SUPPORT PLAN 


The experiment support process can he divided into five distinct phases. 

1. Requirements definition 

2. Requirements implementation 

3. Verification 

4. Production 

5 . Archival 

Each phase represents an important step in the experiment support process, 


2.1 Requirements Definition 

Requirement definition is one of the most important functions of the experir- 
ment support task. In past programs, numerous problems have arisen because of 
incorrect or insufficient definition of PI data requirements. Such problems have 
resulted. ;n delayed delivery of PI data and/or delivery of data that were not 
correct or uhat were of no use to the PI, Many of these problems could have been 
avoided by more complete definition of PI requirements. To avoid repeating past 
mistakes, the initial function of MPB/TRW Task 308, together with the JSC task 
monitor, was to meet with ASTP Pi's and/or their representatives for the purpose 
of obtaining detailed definitions of PI postflight data requirements defined in 
reference 6- These requirements were then documented by means of ASTP Postflight 
Requirements Forms (PRF's) and submitted to the Pi's for review and concurrence 
and ultimate inclusion in the PI contracts. Definition and documentation of 
requirements were the major functions of Task 308 for its first 3 months, July 
through September, 1974, and were continued at a reduced level well into 1975- 


2.2 Requirements Implementation 

Following documentation of the postflight requirements and PI concurrence , 
the process of implementing the requirements began. In many cases, it was not 
necessary to go through a program design phase, which is normally required, since 
many of the parameters requested by Pi’s were available through existing software. 

In particular, the Houston Operations Predict or /Estimator (HOPE) orbit determina- 
tion program and the Apollo Photo Evaluation (APE) photographic analysis program 
included software for computing many of the requested parameters . Certain other 
requested parameters, however, required new software. A new program, known as 
the Apollo-Soyuz Experiment Parameter (ASEP) program, was developed (designed, coded, 
and verified) to compute these new parameters and to act as an executive in the 
merging of all other output data onto a single PI output tape. 


2.3 Verification 

Following implementation of PI requirements, all ASEP outputs were verified, 
both through in-house verification and through a closed-loop PI verification 
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utilizing simulated mission data. In-house verifications consisted of coding and 
logic checks and hand checks of all computer calculations. The PI verification 
data were generated on the basis of an ASTP reference trajectory and reference 
attitude time line (refs. 1 and 2). Each experiment for which data are to be 
provided was active at least once over the intervals chosen for test . Sample 
output data in the form of magnetic tapes, microfilm, and hard copies were shipped 
to each PI for verification of computer compatibility of magnetic tapes and of 
proper computation and display of all requested parameters. 


2 . 4 Production 

Experiment support data production can be divided up into four basic functions 

1. Trajectory data processing 

2. Attitude data processing 

3. Special processing 

4. Output data generation 

Of these functions, the first two can be further divided into the’ preprocessing 
and processing phases. Section 3 describes all the data processing and production 
functions in detail. 


2 . 5 Archival 

All important data required in the postflight experiment support effort will 
be archived for use by the current Pi's, future scientific investigators, and 
other data users. In general, the archival will follow the guidelines enumerated 
below. 

1. All raw data will be archived. 

2. All master PI output products will be archived. 

3. Intermediate data will be archived where appropriate and where no special 
processing is required for the archival. 

4. Data will be archived at the National Space Sciences Data Center (NSSDC). 

All data are being archived for the purpose of giving the PI the option of 
regenerating his support data at some future time. (JSC will have no ASTP support 
capability after January 1, 1976.) Since a probable reason for desiring a regen- 
eration of data is a refinement of mission-specific constants and/or processing 
techniques, the attitude and tracking data are being archived in their raw, uncor- 
rected form. Additionally, all spacecraft ephemeris and refined attitude data 
will be archived. This will allow for a relatively straightforward regeneration 
of experiment -specific parameters, should the PI refine his experiment -related 
constants (for example, instrument misalinements ) . Finally, all PI support prod- 
ucts will be archived to ensure against loss of data and/or bad data tapes. It 
should be noted that the archived PI support tapes will be compatible with the 
Univac 1108 computing system at JSC. Any use of these tapes by Pi's with other 
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computing equipment vill require reformatting of the tapes. Thus, current Pi's 
should take it upon themselves to make certain that they maintain backup copies 
of the computer compatible tapes, to avoid later problems with tape conversions 
For further information on archiving, see reference 8. 

2.6 Schedules 

The bar chart in figure 2-1 depicts the planned time line for the various 
experiment support functions described previously. 
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Figure 2-1.- Experiment support time line. 





3.0 PRODUCTION PLAN 


The purpose of this section is to delineate the experiment data generation 
process: identify the input data requirements, specify the required data pro- 

cessing functions and the associated software, describe the output data products, 
list archival products, and define a tentative data processing schedule. 


3.1 Processing Overview 

The data processing function has been divided into two operations: postflight 

data preprocessing and data processing. Refer to figure 3-1 for a detailed 
illustration of the general processing flow. The specific elements of these 
operations were developed by consideration of the following factors. 

1. Maximize the use of existing software. 

2. Standardize the processing functions and the output products. 

3. Provide for quality control at each level of processing. 

The preprocessing operation consists of processing, editing, interpolating 
(over short time intervals to fill data gaps), and merging the mission input data 
(telemetry data, tracking data, and camera shutter times) in order to produce a 
set of clean, time-ordered, master data files with the use of existing software 
(GREP. GEDIT, HOPE, and LOVE), with data quality checks at each processing step. 

This operation will account for approximately 75 percent of the processing effort 
with respect to labor and calendar time. 

The processing operation consists of the assimilation of mission-specific 
data, the preparation of program control cards, and the generation of experiment 
support data tapes that contain the support parameters in a standard data, record 
format. A nominal set of parameters (state vectors, Sun angles, and gimbal angles) 
will be generated on a one-for-one basis with the input ephemeris/ attitude data, 
whereas special experiment -related parameters will be generated only during the 
specified operation periods associated with each individual experiment. This 
approach is considered the best way to standardize the production process and 
output products. Since the production program, Apollo-Soyuz Experiment Parameter 
(ASEP) , must be developed from "scratch," the design of the program can be modified 
to reflect variations in the processing philosophy. Thus, there is more flexibility 
available for the modification of the production process. It is estimated that 
this operation will require about 25 percent of the labor resources. 

There are five intermediate points in the overall data processing operation 
at which the quality and quantity of the data may be verified. All data listings 
and tapes for the intermediate steps will be retained until the support products 
have been delivered. The retention of these data will aid in the detection and 
isolation of data problems and facilitate the regeneration of data products. 

3.2 Detailed Processing Plan 

The purpose of this section is to provide a more detailed explanation of the 
preprocessing and processing operations. Some of the techniques and procedures 
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Figure 3-1 General processing plan for experiment data support. 













outlined below were originally developed during the Apollo Program for the genera- 
tion of NASA Apollo Trajectory (NAT) tapes and Apollo Photograph Evaluation (APE) 
tapes, which provided ephemeris , attitude, pointing, and field-of-view data for 
numerous lunar orbital experiments. Some of the other techniques and procedures 
that will be used for the ASTP were developed for the reconstruction of 
SKYBET tapes, which provided ephemeris, attitude, and pointing data for Skylab 
orbital experiments. All existing software - GREP, GEDIT, HOPE, and LOVE - was 
developed during the Apollo Program but was utilized on both the Apollo and Skylab 
Programs to provide experiment support data. Thus, the support to be provided 
for ASTP will be based on proven techniques, procedures, and software. 

3.2.1 Preprocessing .- The preprocessing operation is the most complex and 
time consuming of the two operations. The operation can be viewed as four distinct 
functions, as illustrated by figure 3-2. These functions are identified as follows. 

1. Preprocess "raw" gimbal data. 

2. Reconstruct a refined spacecraft trajectory as required. 

•i 

3. Generate trajectory tapes. 

4. Edit, interpolate, and merge trajectory and gimbal data tapes. 

Detailed descriptions are presented below. 

The purpose of function 1 is to strip time and gimbal words from the 
Univac 1108 input tapes, convert the time words into G.m.t. times, eliminate 
duplicate data points, create gimbal data "work" tapes, and list all the gimbal 
data. The listings provide for a check of the quality and quantity of the data, 
as well as the identification of missing data. Once data gaps are identified, 
the missing data may be re-requested or the gaps may be filled by interpolation 
or extrapolation, if possible. Finally, the DELTT cards - HOPE input trajectory 
control cards - may be derived from the data listings in preparation for the 
generation of trajectory tapes. 

Function 2 consists of the preprocessing of the "raw" tracking data tapes 
for input to the orbit determination program. The orbit determination program 
will be used to reconstruct refined trajectories for all the experiments. Special 
techniques, devised by the Mathematical Physics Branch, will be used to reduce 
spacecraft position errors to 200 meters (la), unless there is excessive vehicle 
maneuvering. 

Function 3 consists of the generation of trajectory data tapes for the total 
mission experimental period. There are two distinct modes. The first, a fixed 
frequency mode, generates trajectory data to match exactly the times associated 
with the preprocessed gimbal data. The second mode generates trajectory data 
at the MA-I 36 camera exposure times that have been converted to G.m.t. times. 

The DELTT trajectory control cards for the HOPE program are used for the first 
mode, whereas the SPTME trajectory control cards, derived from the listing of MA-136 
exposure times, are used for the second mode. Refined state vectors will be used to 
generate the ephemerides for all experimental periods * 

Function 4 consists of the editing, interpolating and merging of trajectory 
and gimbal data. For the fixed frequency mode, the gimbal times may be tested 
against the corresponding trajectory times to identify and edit "wild" gimbal 
time points. For short time intervals where gimbal data do not exist, interpolation 
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Figure 3-2 Functional flow for the data preprocessing operation. 











routines will be used to fill in these gaps and such data will be tagged 
accordingly. For the special time mode, the time and gimbal data are linearly 
interpolated in order to match the trajectory times that correspond to the camera 
exposure times. After all editing and interpolation is complete, the data will 
be merged onto master tapes in preparation for the next operation. 

Figure 3-3 depicts the time phasing of these preprocessing functions. Ini- 
tially, functions 1 and 2 may proceed in parallel, as well as some phases of 
function 3. Function 3, however, cannot be completed until all the data associ- 
ated with functions 1 and 2 are generated. In particular, it should be noted 
that function 1 cannot even begin until raw attitude data tapes have been delivered 
by DSAT/IDSD. Function 4 is initiated after the completion of function 3 and 
proceeds until all merged master tapes are assembled. 


Function 

1 

2 

3 

4 


EOM 


Time (months) 


Figure 3-3 Preprocessing operation time phasing. 


3.2.2 Processing .?- The processing operation consists of the computation of 
required ephemeris, attitude, and pointing parameters and the generation of re- 
quired output products - tapes, listings, and microfilm. The operation is 
straightforward once the mission-specific data in the form of cards (see sec. 6.2.1) 
and the master input data tapes are available. The output tape will consist of 
header records of fixed length and data records of fixed length. A nominal set of 
parameters (state vectors. Sun angles, and gimbal angles) will be generated on a 
one-for-one basis with the input ephemeris/attitude data, whereas special experiment 
related parameters will be generated only during the actual operation periods 
associated with each experiment (as specified at end-of-mission (EOM) by the re- 
sponsible Pi’s). If gimbal angle data are not available at a specific time and 
cannot be reconstructed, only ephemeris -related parameters will be computed and 
all attitude and pointing parameters will be flagged with a code word to indicate 
that the gimbal data were missing. The output tape format is specified in 
section 6.2.4. 

With regard to other output products, microfilm will be generated on a one- 
for-one basis with the output data, whereas tab listings will be printed at a 
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frequency controlled "by card input to the program. The format is specified in 
section 6. 2 A. 

3.2.3 Special processing requirements for MA-059 -- To supplement the param- 
eters provided by the Apollo-Soyuz Experiment Parameter (ASEP) program for the 
UV absorption experiment, MA-059, additional parameters will be generated from 
the reconstruction of the relative trajectory with use of the HOPE orbit deter- 
mination program. Special purpose utility programs will be required to preprocess 
the relative tracking data used in the orbit determination process and to merge 
the HOPE parameters with the ASEP parameters. Figure 3-4 displays the required 
processing functions associated with this special support. 


3.3 Computer Compatible Tape Conversions 

The master output data tape generated by the ASEP program, as well as the 
special MA-059 data tape, must be converted to forms compatible with PI computa- 
tional systems. These computer compatible tapes (CCT's) are the tapes that are 
shipped to Pi's. Special purpose utility programs are used for these conversions. 
The tape formats are discussed in appendix C. 

3 . 4 Archival 

It is proposed that the output data tapes (Univac 1108) and associated micro- 
film, as well as the merged data tapes, be archived for possible future use. 

This archival function would have minimal impact with respect to labor and computer 
resources and delivery schedules if this function is specified before the data 
generation is actually initiated. 


3-5 Schedule 

A tentative schedule of data deliveries, based on a current understanding 
of MPAD, IDSD, and principal investigator commitments, is presented in table 3-1, 
where all delivery dates are referenced to the end of mission. 
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Figure 3-4 .- Special-purpose processing functions for MA-059. 
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TABLE 3-1.- DATA PROCESSING SCHEDULE 


Item 

Responsible source 
(organization) 

£L 

Delivery date 

Calandar date 
based on 
nominal mission 

Experiment specific data 

PI/S&AD 

EOM 

7/24/75 

Mission specific data 




As-flown experiment time lines 

PI/S&AD 

EOM 

7/24/75 

Tracking data tapes 

GDSD 

2 weeks 

8/7/75 

Camera shutter times 

IDSD 

4-8 weeks 

8/21-9/18 

Telemetry data tapes 

IDSD 

4-8 weeks 

8/21-9/18 

Mission orbit determination 

MP AD/TRW 

6 weeks 

9/4/75 

Telemetry data rerun request 

MPAD/TRW 

9 weeks 

9/25/75 . 

Merged gimbal & trajectory tapes 

MP AD/TRW 

10 weeks 

10/1/75 

Experiment support data tapes 

MPAD/TRW 

14 weeks 

11/1/75 

Experiment support data rerun request 

PI 

18 weeks 

11/29/75 

Final experiment support data tapes 

MPAD/TRW 

23 weeks 

12/31/75 

Archival 

MPAD 

34 weeks 

3/15/76 

Final report 

MPAD 

34 weeks 

3/15/76 


a All dates referenced to the end of mission. 



4.0 EXPERIMENT SUPPORT PARAMETER FORMULATIONS 


The purpose of this section is to specify the equations for all the experi- 
ment support parameters that will he provided by JSC/TRW Task 308. The 
following subsections define the standard nomenclature for all formulations , define 
all required computational coordinate systems and associated transformations, and 
specify the equations for all the experiment support parameters. 


4 . 1 Nomenclature 

This section defines the nomenclature to be used for the specification of 
experiment parameter formulations. The units associated with the parameters are 
kilometers, seconds, and degrees unless otherwise specified. The coordinate 
system numeric designators are listed below. The coordinate systems are defined 
in section 4.2. 


Numeric 

Definition identi fier 

Mean of 1950.0 1 

Ecliptic mean of 1950.0 2 

True of date 3 

Geographic (inertial) 4 

IMU (nominal) 5 

IMU ( actual ) 6 

Navigation base 7 

Instrument base (nominal) 8 

Instrument base (actual) 9 

UEN - up, east, north 10 

Geomagnetic 11 

Geographic (rotating) 12 


3 x 3 transformation matrix from the nominal instru- 
ment base (8) to the actual instrument base (9)* 

Semimajor axis of the Earth’s reference ellipsoid. 

ASTP mission elapsed time from range zero. 

Instrument field-of-view designator. 
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A 

a 

AET 

Ai 




AL 

ALFi 

ALPHA 

ALPH10, BETAIO , 
PHIIO 

ALT4 

ALTR 

AZi 

B 

b 

Bi 

BMAG 

BTAi 

C 

CDUX, CDUY, CDUZ 

Ci 

CET 

D 

Di 


sSKawsaaEi- 


Surface arc length between nadir and the principal 
point. 

Right ascension of the vehicle measured from the 
x-axis of the i^h coordinate system. 

Angle between the camera z-axis and the projection 
of the local vertical (principal point) into the 
plane formed by the vector from the principal point 
to the Sun and the camera z-axis. 

Three orientation angles that relate the spacecraft 
body axes to the UEN coordinate system. 

Geodetic altitude of the vehicle. 

Time rate of change of the altitude (inertial) of 
the spacecraft with respect to the principal point. 

Azimuth of the velocity vector in the i th coordinate 
system measured clockwise from the local north to 
the projection of the velocity vector in a plane 
normal to the local vertical. 

3 x 3 transformation matrix from the geographical 
inertial coordinate system (4) to the geomagnetic 
coordinate system (ll). 

Semiminor axis of the Earth's reference ellipsoid. 

Instrument field-of-view designator. 

Intensity of the Earth's magnetic field (gauss). 

Vertical flightpath angle measured from the local 
vertical to the velocity vector in the i^ h coordinate 
system. 

3 x 3 transformation matrix fr^-n the geographical 
inertial coordinate system (4) to the UEN coordinate 
system ( 10 ). 

Outer, inner, and middle gimbal angles. 

Instrument field-of-view designator. 

Onboard clock elapsed time for the CSM central 
timing equipment. 

3x3 transformation matrix from the nominal 1MU 
coordinate system (5) to the actual IMU coordinate 
system ( 6 ). 

Instrument field-of-view designator. 
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DLTi 


DT 

E 

ECCi 

EMISS 

ESLOG 

ESLOSS 

F 

f 

F3RHA1 , 
F3RHB1 , 
F3RHC1 , 
F3RHD1, 
F8RHA1, 
F8RHB1 , 
F8RHC1 , 


Declination of the vehicle in the 1 th coordinate 
system. 

Time interval between adjacent photographs. 

3x3 transformation matrix from the mean of 1950.0 
coordinate system (l) to the ecliptic mean of 1950.0 
coordinate system (2) (also denoted E.,). 

Eccentricity of the orbit in the i^* 1 coordinate 
system. 

Angle between the camera optical axis and the local 
vertical at the principal poirit. 

Angle between the Earth-Sun vector and the SIM Bay 
sector 1 (MA-048, MA 083, MA-088) line-of-sight . 

Supplement of ESLOS. 

3 x 3 transformation matrix from the navigation 
base coordinate system ( 7 ) to the nominal instrument 
base coordinate system (8). 

Earth's flattening (unitless). 

F3DCA1 Right ascension and declination associated with the 

"A quadrant" of the MA-083 field-of-view in the 
mean of 1950*0 coordinate system (l). 

F3DCB1 Right ascension and declination associated with the 

"B quadrant" of the MA-083 field-of-view in the 
mean of 1950.0 coordinate system (l). 

F3DCC1 Right ascension and declination associated with the 

"C quadrant" of the MA-083 field-of-view in the 
mean of 1950.0 coordinate system (l). 

F3DCD1 Right ascension and declination associated with the 

"D quadrant" of the MA-083 field-of-view in the 
mean of 1950.0 coordinate system (l). 

F8DCA1 Right ascension and declination associated with the 

"A quadrant" of the MA-088 field-of-view in the 
mean of 1950.0 coordinate system (l). 

F8DCB1 Right ascension and declination associated with the 

"B quadrant" of the MA-088 field-of-view in the 
mean of 1950.0 coordinate system (l). 

F8DCC1 Right ascension and declination associated with the 

"C quadrant" of the MA-088 field-of-view in the 
mean of 1950.0 coordinate system (l). 



F8RHD1, 

FL 

FOV 

G 

GC 

GET 

GHA 

GIMB 

GMT 

HEAD 

HV 

INCi 

L1AD1 

L1AR1 

L1AZ 

L1DEC1 

L1DEC2 

L1EL 


F8DCD1 Right ascension and declination associated with the 

"D quadrant" of- the MA-088 f ield-of-view in the 
mean of 1950.0 coordinate system (l). 

Focal length of the camera lens. 

Field-of-view. 

3x3 transformation matrix from the actual IMU 
coordinate system (6) to the navigation "base 
coordinate system (j). 

3 * 3 transformation matrix from the geographic 
coordinate system (4) to camera axes coordinate 
system. 

Ground elapsed time . 

Angle "between the vernal equinox and the Greenwich 
meridian. 

Gimbal data status flag. 

Greenwich mean time (UTC). 

Angle, measured positive clockwise in the local- 
horizontal plane, from local north to the projection 
of the camera x-axis into the local-horizontal plane. 

Component of the spacecraft ' s velocity vector 
(inertial) that is colinear with the local-horizontal 
plane with respect to the principal point . 

Inclination of the orbit plane with respect to the 
x-y plane in the i^h coordinate system. 

Alinement corrected angle L1DEC1. 

Alinement corrected angle L1RHA1. 

Azimuth of the SIM Bay sector 1 (MA-048, MA-083, 

MA-088) line-of-sight in the UEN coordinate 
system (10) . 

Declination of the SIM Bay sector 1 (MA!-048, MA-083, 
MA-088) line-of-sight in the mean of 1950.0 co- 
ordinate system (l). 

Declination of the SIM Bay sector 1 (MA-048, MA-083, 
MA-088) line-of-sight in the ecliptic mean of 
1950.0 coordinate system (2). 

Elevation of the SIM Bay sector 1 (MA-048, MA-083, 
MA-088) line-of-sight in the UEN coordinate system (10) 


L1RHA1 

L1RHA2 

LAT4 

LATA, LOWA 

LATB, LONB 

LATC, LONC 

LATD, LOUD 
LATP, LOUP 
LATS, LOWS 

LH 

LMAG 
LON4 
LONS 2 

LOS 

LOSX, LOSY, 
ni 

NODi 


LOSZ 


Right ascension of the SIM Bay sector 1 (MA--048, 
MA-083, MA.-O 88 )- line-of-sight in the mean of 
1950.0 coordinate system (l). 

Right ascension of the SIM Bay sector 1 (MA-048, 
MA.-O 83 , MA- 088 ) line-of-sight in the ecliptic mean 
of 1950.0 coordinate system ( 2 ). 


Geodetic latitude of the vehicle. 

Geodetic latitude and longitude associated with the 
"A corner” of the MA-136 field-of-view. 


Geodetic latitude and longitude associated with the 
”B corner” of the MA-136 field-of-view. 


Geodetic latitude and longitude associated with the 
"C corner” of the MA-136 field-of-view. 

Geodetic latitude and longitude associated with the 
”D corner" of the MA-136 field-of-view. 

Geodetic latitude and longitude of the principal 
point . 

Latitude and longitude of the intersection of the 
vector from the Earth’s center to the Sun with the 
Earth's ellipsoid. 

3x3 transformation matrix from the geographic 
inertial coordinate system (4) to UEN coordinate 
system (10) . 

Radius of the magnetic field ”L shell” (Earth radii) 

Geographic longitude of the vehicle. 

Longitude of the Sun in the ecliptic mean of 1950.0 
coordinate system (2). 

Line-of-sight . 

Direction cosines for a vector from the vehicle to 
the principal point in the geographic inertial 
coordinate system (4). 

Subscript that designates the instrument line-of- 
sight or the camera optical axis in the i™ co- 
ordinate system. 

Longitude of the ascending node, measured in the 
x-y plane from the x-axis in the ith coordinate 
system. 
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OMGi 


OPFLAG 

OVR 

PET 

PHASE 

PHI, KAPPA, 
OMEGA 


pi 


PP 


psi 


R 


r 

REV 

RF 

Ri 

SAZIO 

SAZP 


Argument of perifocus, measured in the orbit plane 
from the ascending node to the point of closest 
approach in the 1 th coordinate system. 

Flag that designates which experiments are operating. 

Ratio of the amount of overlap of the photograph 
frame dimension along the direction of the flightpath. 

Phase elapsed time. 

Angle between the camera optical axis and the. vector 
from the principal point to the Sun. 

Angles that rotate the camera axes coordinate system 
into the local-horizontal coordinate system, where 

<f> = primary right-handed rotation about the camera 
y-axis . 

w = secondary right-handed rotation about the inter- 
mediate x-axis. 

< = final right-handed rotation about the local- 
vertical axis. 

Subscript that denotes a vector from the Earth's 
center to principal point in the ith coordinate 
system. 

Principal point - point of intersection of the 
camera optical axis with the Earth's ellipsoid. 

Subscript that denotes a vector from the principal 
point to the Sun in the i^* 1 coordinate system. 

3*3 transformation matrix from the mean of 1950.0 
coordinate system (l) to the nominal IMU coordinate 
system ( 5) • 

Vehicle position vector. 

Vehicle rev. number. 

3*3 transformation matrix that relates the stable 
member axis or nominal IMU ( 5 ) system to the mean 
of 1950.0 coordinate system. 

Magnitude of the position vector. 

Azimuth of the Sun vector in the UEN coordinate 
system (10) . 

Sun azimuth at the principal point - angle, measured 
positive clockwise, from local north to the projection 
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SCSA 

SEL10 

SELP 


SF 

si 

SMAi 

SR 

SUNF 

SUNH, SUNM, 
SUNS 

SWING 

T 

TAi 

THETA 

TILT 

TILTAZ 

v 


of the vector from the principal point to the 
Sun into the local-horizontal plane. 

Angle between the Earth- centered vehicle position 
vector and the Earth-Sun vector. 

Elevation of the Sun vector in the UEN coordinate 
system (10). 

Sun elevation at the principal point - angle, 
measured counterclockwise, between the vector from 
the principal point to the Sun and the local- 
horizontal plane at the principal point. 

Constant that relates film dimensions to surface 
dimensions. 

Subscript that denotes a vector from the Earth to 
the Sun in the ith coordinate system. 

Semimajor axis of the orbit in the i*^ coordinate 
system. 

Slant range from the vehicle to the principal point. 
Sunrise/sunset designator. 

Sunrise/sunset time. 


Angle, measured counterclockwise, between the camera 
y-axis and the projection of the local vertical 
into the camera x-y plane. 


3 x 3 transformation matrix from the mean of 1950.0 
coordinate system (l) to the true of date coordinate 
system (3) (also denoted T^ ) . 


"t h 

True anomaly in the i coordinate system, measured 
in the orbit plane from the point of closest approach 
to the vehicle's position in the orbit. 


Angle between the vehicle x-axis and the CSM velocity 
vector. 

Angle between the camera optical axis and the local 
vertical. 

Angle, measured clockwise, between local north and 
the projection of the camera optical axis into the 
local-horizontal plane. 

Vehicle velocity vector. 
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Vi 

vi 

VDEC1 

VDEC2 

VLOS 

VRHA1 

VRHA2 

W 

x 

Xi, Yi, Zi 
XANG 

XDi, YDi , ZDi 


Vehicle speed. 

Subscript that denotes a vector from the center 

"fcll 

of the Earth to the vehicle in the i coordinate 
system. 

Declination of the vehicle velocity in the mean 
of 1950.0 coordinate system (l). 

Declination of the vehicle velocity in the ecliptic 
mean of 1950.0 coordinate system ( 2 ). 

Angle between the vehicle velocity vector and the 
instrument line-of-sight . 

Right ascension of the vehicle velocity in the 
mean of 1950.0 coordinate system (l). 

Right ascension of the vehicle velocity in the 
ecliptic mean of 1950.0 coordinate system (2). 

3x3 transformation matrix from the true of date 
coordinate system ( 3 ) to the geomagnetic inertial 
coordinate system (4). 

Denotes an instrument pointing vector. 

■fch 

Vehicle position in the i Cartesian coordinate 
system. 

x-angle deviation of the vehicle-Sun vector with 
respect to the instrument line-of-sight. 

"til 

Vehicle velocity in the i Cartesian coordinate 
system. 


XDSi ,YDSi , Velocity of the Sun with respect to the Earth in the 

ZDSi mean of 1950.0 coordinate system (l). 


XSi , YSi , ZSi 
XTIDT 
YANG 
YTILT 
« 1 » <* 2 , « 3 


Position of the Sun with respect to the Earth in the 
mean of 1950.0 coordinate system (l). 

Angle from the local-horizontal plane to the camera 
y-axis (lateral tilt). 

y-angle deviation of the vehicle-Sun vector with 
respect to the instrument line-of-sight . 

Angle from the local-horizontal plane bo the camera 
X-axis (longitudinal tilt). 

Inner, middle, and outer gimbal angles. 
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4.2 Definition of Coordinate Systems 

A total of 12 coordinate systems is involved in the computation of the required 
experiment support parameters. The purpose of this section is to define each 
required coordinate system. 

4.2.1 Mean of 1950.0 .“ The x-axis of this inertial Cartesian coordinate system 
lies along the line of intersection of the mean celestial equator with the ecliptic 
plane at 1950.0 and points toward the vernal equinox. The z-axis is orthogonal to 
the mean celestial equator at the reference time, and the y-axis completes the 
right-handed system. 

4.2.2 Ecliptic of 1950»0 »- The x-axis for this inertial Cartesian coordinate 
system lies along the line of intersection of the mean celestial equator with the 
ecliptic plane at 1950.0 and points toward the vernal equinox. The z-axis is 
orthogonal to the ecliptic plane at the reference time, and the y-axis completes 
the right-handed system. 

4.2.3 True of date .- The x-axis of this inertial Cartesian coordinate system 
lies along the intersection of the true celestial equator of date and the ecliptic 
plane and points towards the vernal equinox of date. The z-axis is orthogonal to 
the true celestial equator of "date, and the y-axis completes the right-handed 
system. 


4.2.4 Geographic inertial .- The x-axis of this Cartesian coordinate system 
lies in the celestial equator of date through the Greenwich meridian of date. The 
z-axis is orthogonal to the celestial equator of date, and the y-axis completes the 
right-handed system. 

4.2.5 IMU ( nominal ) . - The CSM IMU Cartesian coordinate system is alined 

in inertial space to the mean of 1950.0 system through the Apollo REFSMMAT, which 
relates the stable member (SM) axis system to the mean of 1950.0 coordinate system. 
The nominal system is assumed to he error free; that is, no misalinements , drifts, 
et cetera. 


4.2.6 

nominal IMU 


IMU ( actual ) . - The actual IMU coordinate system is related to the 
system by means of platform drifts. 


4.2.7 Navigation base .- The navigation base coordinate system is a Cartesian 
system that corresponds to the vehicle-fixed x , y^, and z^ coordinate system. 

4.2.8 Ins trument base (nominal) .- For each instrument, as required, a 
Cartesian coordinate system is defined such that the z-axis is colinear with the 
nominal (perfectly alined) instrument LOS and orthogonal to the x-y plane. The 
specification of the x, y axes will be instrument dependent . 


4.2.9 Instrument base (actual) .- This Cartesian coordinate system is related 
to the nominal instrument base coordinate system by means of three-axes rotation 
matrix of misalinement angles. 


4.2.10 Up, east, north (UEN) .- The x-axis of this local-horizontal coordinate 
system is colinear with the position vector of the vehicle. The z-axis is directed 
towards local north in the geographic inertial system, and the y-axis (directed 
towards local east) completes the right-handed system. 
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4.2.11 Geomagnetic The z-axis of this Cartesian coordinate system is 
parallel to the extended axis of a "best fit” of a centered dipole to the Earth's 
actual magnetic field, the y-axis lies in the celestial equator of date, and the 
x-axis completes the right-handed system. The location of the intersection of the 
extended axis of the centered dipole with the Earth's surface is +69.5° west longi- 
tude and 78*5° north latitude. 

4.2.12 Geographic rotating .- This Cartesian coordinate system is identical to 
the geographic inertial coordinate system except that it is rotating about the 
z-axis with angular speed ^ e , where w e is the rotation rate of the Earth. 

4.2.13 Local vertical .- The local-vertical coordinate system, although not 
used directly in experiment support computations, is useful and often referenced 
for analysis purposes. This Cartesian system is Earth- centered . Its x-axis 
(radial axis) is directed radially outward along the vehicle position vector. The 
z-axis (crosstrack axis) is in the direction of r x v, where r is the vehicle 
position vector and v is the velocity vector. The y-axis (downtrack axis) com- 
pletes the right-handed system. 

4.3 Coordinate System Transformations 

To compute the required experiment support parameters, it will.be necessary, 
in general, to transform vectors from any one of 12 possible coordinate systems, 
as defined in section 4.2, to any one of the remaining coordinate systems. The 
definition of these transformations can be simplified by the visualization of three 
"strings” of transformations that originate from the base system (mean of 1950.0). 
These strings of transformations include vehicle related transformations, true 
Earth equatorial related transformations, and the ecliptic mean of 1950.0 transfor- 
mation, as illustrated in figure 4-1. Thus, any required transformation can be 
computed on the basis of the proper selection of matrix products as indicated in 
the figure. 

The objective of this section is to specify the formulation for those trans- 
formations that must be computed or to identify the source of those transformations 
that will be defined by input (tape or card). 

4.3.1 Mean of 1950.0 to nominal IMU .- The transformation R is specified by 
the REFSMMT matrix through card input . 

4.3.2 Nominal IMU to actual IMU .- The transformation D is computed as follows, 
on the basis of the drifts that are defined in a right-handed sense about the x, y, 
and z axes . 


D = <[l] - ^ sin (d At) [H] + [l - cos(d At)] [H] 2 
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Figure 4-1." Coordinate system transformation "tree. 









where 


[i] = identity matrix, 3 * 3 

d: (d , d , d ) T = drift about the IMU axes 

x* y’ z 

d = (d-d) 1/2 

At = time elapsed since alinement of the IMU 


[H] = 


0 -d d 
z y 


d 0 -d 
z x 


-d d 0 

y x 


The values for d , d , and d are obtained through card input, 
x’ y’ z 

4.3.3 Actual IMU to navigation base .- The transformation G is a function of 
the inner (a^ ) , middle (a 2 ), and outer (a^) gimbal angles, as follows. 

G = [g iJ ] ; i,j =.l, 2, and 3 

g = cos a. cos a 0 
11 12 

g 12 = Sin a 2 

g 13 = -sin cos a 2 

g 21 = -cos sin a 2 cos + sin sin 


"22 

5 23 

"31 


= cos 

CVJ 

a 

cos 

a 3 







= sin 

a l 

sin 

a 2 

cos 

a 3 

+ 

cos 

sin 

a 3 

= cos 

“l 

sin 

a 2 

sin 

“3 

+ 

sin 

cos 

a 3 


g 32 = -cos « 2 sin a 3 


"33 = “ s: *- n a l sin a 2 sin a 3 + cos “i cos a 3 


The gimbal angles are available through tape input. 
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4.3.4 Navigation base to nominal instrument base .- The transformation F is 
instrument dependent; it is defined in section' 4.4 as required. 

4.3.5 Nominal, instrument base to actual instrument base .- The transformation 
A is instrument dependent; it is defined in section 4.4 as required. 

4.3.6 Mean of 1950.0 to ecliptic mean of 1950. 0 .- The transformation E is 
available through the tape input as computed by the HOPE program (see app. A). 

4.3.7 Mean of 1950.0 to true of date .- The transformation T is available 
through the tape input as computed by the HOPE program (see app. A). 

4.3.8 True of date to geographic inertial . - The transformation W is computed 
as follows . 


W = (w^ ) ; i,j » lj 2, and 3 

where 

V 11 = W 22 = COS ° 
w 2i = -v = -sin a 

W 31 = V 13 = W 32 = W 23 = ° 

W 33 = 1 

a = a + o) At 
go e 

= right ascension of Greenwich, midnight day of epoch 
oj = Earth's rotation speed 

At = elapsed time since midnight day of epoch 


The value for a 

go 


is available by tape input. 


4.3.9 Geographic inertial to up, east, north (UEN) .- The transformation C is 
computed as follows. 


C = "(c ). ; i,j = l, 2 and 3 
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where 


c^ = cos <j> cos X 


c g = cos $ sin X 


C 12 = Sin 


c = -sin X 


c 22 ~ cos X 


°23 = 0 


c^ ~ -sin <f> cos X 


C32 = -sin <p sin X 


C — cos 

33 


and where <f> and X are the computed geodetic latitude and longitude, respectively. 


*+•3.10 Geographic inertial to geomagnetic .- The transformation B is computed 
as follows. 


B = (by) 1 U = 1 » 2 » and 3 


where 


b^ = cos y cos X 


t> 12 = cos y sin X 


b 13 = -sin Y 


h 21 = -sin X 
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b 22 - cos X 



b = sin y c os X 
bg 2 = sin Y sin * 

*33 = 008 * 

and where X = -69.5° and y = 11.5°, the longitude and colatitude, respectively, of 
the Earth’s magnetic North Pole. 

4.3.11 Geographic inertial to geographic rotating .- The transformation of a 
state vector from the geographic inertial coordinate system to the geographic ro- 
tating coordinate system is defined as follows. 

? 12 = ? 4 

*12 " *4 - “ * ? 4 

5 = ( 0 , 0, a ) T 

= rotational speed of the Earth 
4.4 Experiment Support Parameters 

The objective of this section is to specify in detail the formulation for. each 
experiment support parameter to be provided to ASl'P principal investigators. There 
are two basic sets of parameters: the nominal ephemeris and pointing parameters 

and the experiment-specific (attitude, pointing, and field-of-view) parameters. 
Section 4.4.1 specifies the nominal parameters, whereas sections 4. 4. 2.1 through 
4. 4.2. 8 specify the experiment-related parameters. 

4.4.1 Nominal support parameters .- The set of nominal support parameters will 
be divided into four categories, as follows. 

1. Time-related parameters. 

2. Vehicle state vectors . 

3. Vehicle-attitude-related parameters, 

4. Miscellaneous support parameters. 
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These categories will be discussed in the above-mentioned order in the follow- 
ing subsections. 

4. 4. 1.1 Time-related parameters: The six sets of time-related parameters are 

specified as follows. 

1. G.m.t. - Greenwich mean time (UTC). 

2. g.e.t. - Ground elapsed time. The base time for g.e.t. may be updated 
during the mission. 

3. AET - ASTP mission elapsed time from range zero. 

4. CET - Onboard clock elapsed time. All onboard computer parameters are 
tagged by CET. 

5. PET - Phase elapsed time, which is used to designate special experiment 
phases . 

6. OPFLAL - Flag to designate which experiments are operating. 

All parameters will be displayed in terms of hours, minutes, and seconds from a 
base time. 

The OPFLAG is a seven digit number used to designate experiment operation 
times. Each digit refers to a specific experiment, as follows. 


Digit ft 

Experiment ft 

Experiment 

1 

1 

MA-007 

2 

2 

MA-048 

3 

3 

MA-059 

4 

4 

MA-083 

5 

5 

MA-088 

6 

6 

MA-106 and/or MA-107 

7 

7 

MA-136 

When an experiment is inactive, its digit 
experiment-specific parameters will be computed 

in the OPFLAG will be zero, 
for that experiment. If the 


periment is active, its digit in the OPFLAG will equal its experiment number as 
shown above. Thus, if MA-007 and MA-106 are active while all other experiments are 
off, the OPFLAG will be IOOOO 6 O.E+OO. Similarly, if MA-048, MA-083, and MA -088 are 
active, the OPFLAG will be 0204500 .E+00. 
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4. 4. 1.2 Vehicle state vectors: The vehicle state vectors will he displayed 

in four reference coordinate systems and three' different element sets, as designated 
below. 

Coordinate systems Element sets 

Mean of 1950-0 Cartesian 


Geographic rotating 

True of date Spherical polar 

Geographic inertial Keplerian 

The coordinate systems have been specified in section 3.1*. 2; the element sets are 
defined below. 

Cartesian elements 

• * • 

Th^ Cartesian elements x, y, z, x, y, and z are the components of the state 
vector r, v of the vehicle at time t in the requested output coordinate system, 
with r and v represented by the column vectors . 

(x, y, zj 

/* • ■ \T 

(x, y, z} 

The output elements are denoted X, Y, Z, XD, YD, and ZD. 

Spherical elements 

The spherical elements ALF, DLT, BTA, AZ, R, and V are determined from the 
state vector r, v at time t as follows. 

til 

ALF is the right ascension of the vehicle measured from the x-axis in the i 
coordinate system. 

ALF = tan” 1 (fy ; 0° < ALF < 360° 

DLT is the declination of the vehicle. 

DLT = tan 

BTA is the vertical flightpath angle measured from the local vertical to the 
velocity vector v. 


O „ . I„ ; -90° < DLT < 90° 
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BTA = cos 


1 (^-) 


0° < BTA < l8d° 


where r - jr] and v = jv| 


AZ is the azimuth of the velocity vector measured clockwise from the projec- 
tion of the reference z-axie to the projection of the velocity vector in a 
plane normal to the local vertical. 


AZ = tan 


i 0° <_ AZ < 360 c 

\rz - zrl 


i i 

4, 


where 


• _ r*v 

r 


R is the magnitude of the position vector r. 


/ 2 ^ 2 . 2VL/2 

= (x + y + z J 
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V is the magnitude of the velocity vector v. 

/ *2 *2 *2X1/2 

V = I x + y + z I ' 

Keplerian elements 

The Keplerian elements SMA, ECC, INC, NOD, OMG, and TA are determined from the 
state vector r, v of the vehicle at the time t as follows. 

SMA is the semimajor axis of the orbit. 

a > 0, elliptical orbit 
a < 0, hyperbolic orbit 
a = 0, parabolic orbit 


SMA = 


yr 


2y - rv 


where 


y = gravitational parameter of the Earth 


ECC is the eccentricity of the orbit. 


ECC = (1 - p/a) 1//2 0 <_ ECC < 1 , elliptical 

ECC < 1 , hyperbolic orbit 
ECC = 1 , parabolic orbit 

where 

22 /- -v2 

p = r v - (r»v) 

y 

INC is the inclination of the orbit plane with respect to the x-y plane in the 
coordinate system,. 

INC ~ tan 


_(yz - zy) + (zx - xz)' 

(sty - yx) 


1/2 


INC < 180° 
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NOD is the longitude of the ascending node measured in the x-y plane from the 
,th 


x-axis in the i coordinate system. 

‘‘ [-aV-tz ' j 360 ° 


NOD = tan 


OMG is the argument of perifocus measured in the orbit plane from the 
ascending node to the point of closest approach. 


OMG = tan" 1 


(up) 1 / 2 z 


y(yz - zy) - x(zx - xz) 


- v ; 0° < OMG < 360° 


where 


v = true anomaly 


TA is the true anomaly measured in the orbit plane from the point of closest 
approac . to the vehicle's position in the orbit. 


TA = tan 1 


(r-v) (up ) 1|/2 
(up - up) 


; 0° <_ TA < 360° 


4. 4. 1.3 Vehicle-attitude-related parameters: The vehicle-attitude-related 

parameters are specified as follows. 

1. RF(REFSMMAT) - A 3 x • 3 matrix that relates the orientation of the IMU 
stable member axes system to the mean of 1950 coordinate system. This matrix 
is available by means of card input. 

2. CDUX, CDUY, CDUZ - The inner, middle, and outer gimbal angles that relate 
the vehicle navigation base coordinate system to the IMU stable member axes 
system. These angles are available by means of tape input. 

3 . GIMB - Gimbal data quality flag that denotes the status of the parameters 
CDUX, CDUY, and CDUZ. 


GIMB 


0 Good data 

1 Interpolated data 

2 No data 

3 Bad data 

4 Premission attitude reference data (MA-048 only) 


This flag is available through tape input. When no gimbal data are available or 
the data are bad, the value of all attitude, pointing, and FOV parameters is 
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0.77777770+7* For the MA-136 experiment, a value of 0.888Q8888+8 for the geodetic 
location of an FOV corner point (s) or the principle point indicates that the 
associated line-of-sight is "above the horizon." If this is the case, then any 
other parameter that depends on the principle point or the FOV corner point (s) 
will be set equal to 0.88888888+8. 


ll.U.l.i* Miscellaneous support parameters: The miscellaneous parameters 

that are to be computed are described below. 

Sun angles 

The Sun angles are defined to be the azimuth and elevation of the vector from 
the subvehicle point to the Sun (r'^g) with respect to the UEN coordinate system. 


r ’l0 ■ c(mr sl - 


SEL10 = tan"' x' lQ /(y' 1Q 2 + 10 2 '> 1,S 


-00° < SEL10 < 90 


SAZ10 = tan" 1 ( Y'^q/ 2 


1 (*’io /, 'io) : 00 i 


SAZ10 < 360° 


r’> - subvehicle vector = (r cos X cos <}> , r sin X cos <J> , r sin <jt)‘ 


r = r i, - ALT 
a vt 


Subvehicle point 

The subvehicle location is defined as follows. 


LAT is the geodetic latitude of the vehicle defined by_ the projection of the 
vehicle’s position perpendicular to the reference ellipsoid. 


LAT = tan 


-1 


% 


(*g + y|) 1/2 ( l - f ) 2 


i -90° < LAT < 90° 


LON is the longitude of the vehicle. 


LON = tan Ky^/x^) ; 0° <_ LON < 3^0° 


ALT is the altitude above the reference ellipsoid for the Earth. 


3 >4 


a(l - f) 


i 



where 

r = magnitude of the position vector of the vehicle 
f = 1 - b/a, flattening of the Earth 
a = semimajor axis of the Earth's reference ellipsoid 
b = semiminor axis of the Earth's reference ellipsoid 

x 4 , y^, z 4 = geographic inertial Cartesian position coordinates of the vehicle 
Greenwich hour angle 

The angle between the vernal equinox and the Greenwich meridian in the true 
of date coordinates at midnight of base day is available through tape input. 


Sun position and velocity 

The Cartesian coordinates of the Sun with respect to the Earth (XSl, YS1, ZS1, 
XDS1, YDS1, and ZDSl) in the mean of 1950.0 are available through tape input. 


Coordinate transformations 

Two required coordinate transformations that are computed by the HOPE program 
and that are available through tape input are as follows. 

T - Mean of 1950.0 to true of date transformation (3x3 matrix). 

E - Mean of 1950.0 to ecliptic mean of 1950.0 transformation (3x3 matrix). 



Rev ■ number 

The revolution number of the vehicle is computed with reference to the launch 
longitude . 


REV = REVI + n 


where REVI (integer) is the initial rev. number input by card and n is the number 
of times the vehicle has passed the launch longitude since REVI (not to be con- 
fused with orbit number, which is based on ascending node crossings). 


** 

i ■ 1 ) 
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4.4.2 Experiment-specific support parameters .- The formulations for the 
experiment-specific parameters are discussed in the following order. 


1 . 

MA-007 

- Stratospheric Aerosol Measurement 

2. 

MA-048 

- Soft X-Ray 

3. 

MA-059 

- UV Absorption 

4. 

MA-083 

- Extreme-UV Survey 

5. 

MA-088 

- Helium Glow 

6. 

MA-106 

- Light Flash 

7. 

MA-107 

- Biostack 

8. 

MA-136 

- Earth Observations 


For a definition of each experiment-specific support parameter, refer to section 4.1. 

The matrix F, the navigation "base coordinate system to the nominal instrument 
base coordinate system transformation, is required for three experiments - MA-007) 
MA-083, and MA-136. The formulations for these experiment-specific matrices will 
he specified in the appropriate section. For all other experiments, F is assumed 
to he the identity matrix. It is also assumed that the matrix A, the nominal to 
actual instrument base transformation, is the identity matrix for all experiments 
unless other data are provided by the PI. 

4. 4. 2.1 Support parameters for the stratospheric aerosol measurement experi- 
ment (MA-007): The additional support parameters for the stratospheric aerosol 

measurement experiment are defined below. Since the MA-007 experiment requires 
line-of-sight angular deviations, it is necessary to define the navigation base to 
nominal instrument base coordinate transformation F. The transformation F for 
MA-007 is defined as follows. , 


x 



where 

•v rp 

z — (cos 8, sin 0 sin <f > , -sin 0 cos 4)' 
y = zx (l, 0, 0) T 

y/|y| 

y * z 
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The nominal values for 0 and are 57-5° and 58.9° > respectively. Refer to 
figure 4-2 for an illustration of the nominal instrument base coordinate system. 


Sunrise/sunset time (SUNK, SUNM, and SUNS) 

The vehicle, Earth, and Sun geometry associated with the sunrise /sunset 

calculation is depicted, on an exaggerated scale, in figure 4-3. The procedure 

is to test for the approximate conditions of sunrise/sunset (within 0.25 second) 

with the, use of a mean radius of the Earth (r ). The final calculation is based 

em 

on the proper geodetic Earth radius, r . The effects of refraction will not be 
considered . 

From the figure, the angles 8 and S can be approximated as follows. 

r 

B m = c ° s ' 1(r « / lv L i> 

4 V = tan _1 [(P r - r em > /r j 

S = radius of the Sun 
r 

r = mean radius of the Earth 
em 


r = mean distance to the Sun 
s 


where 8 is accurate to within 0.0164 degree and AS ’ is accurate to within 
m r 

0.25 x 10™ 3 degree. The condition for sunrise or sunset is 


0=8+ AS =8 + AS ' 

r m r 


Thus , when 


9 i = * /2 - oos " 1 ! - r r,l' r vl / ' r vl' ) i- 0 i 0 i+l 


where 


r sl = r S-d r Bl 


i = designates time t^ 
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Figure 4-2.- MA-007 coordinate system for the right side window. 
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Figure 4-3 Geometry for the sunrise and sunset computation. 
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then the sunset time may he calculated as follows . 

*« ■ *i + i - ( e i + i - 8 ) >K 

= ( e i + i - - *i) 

0 - B + AS 

r 

6 = cos’* [r e /|F vl (t i+1 )]| 

AS r = tan'' [(S r - r e )/|r sl (t 1+1 )|] 

where r is computed as presented below, 
s 

In a similar manner, the sunrise time may be computed as 



where 


e.^ n ifl is. 

i+i i 


The geodetic radius r^ may be computed in the following manner 


r sl * r sl / l r J 

r vl r vl / I r vl ^ 


n = r . x r / Ir.xr. 
s si vl 1 si vl 1 


n = n x r . 
r s si 


S * = 
r 


-S if r *n < 0 
r vl r 

S if r *n > 0 
r- vl r 


r , = r ., + S *n - r _ 
svl si r r vl 


ko 



r svl r svl/^ r svlJ 


r .=r/ir|=n x r , 
el e ' e 1 s svl 


IK (S *)] 
r 


r e4 = r el 


If r is given, the geodetic radius may be computed directly. 


where 


a (l - f) 

_ __e 

[1 - f{2 - f) (x 2 + y 2 )]' /2 


a e = Earth ’ s aa “atorlal radius 
f = Earth's flattening 

r e 4 = < x > y» z ) 


Funrise/sunset flag (SUNE) 


A flag to denote whether the listed time is associated with the vehicle 
sunrise time or the vehicle sunset time is included in the narameter list. The 
flag will be defined as follows . 


0 Not being used 
SUNF = | 1 Sunrise 

2 Sunset 


L ine-of-sight angular deviations (XANG and YANG) 

The coninutations for the required angles are as follows. 


r vsl ' r sl “ r vl 
r vsl r vsl^ | r vsl| 
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= AFGDR r 


vs9 

XANG 

YANG = 


vsl 


* tan' 1 (WK s9 !) i -90° £XANG < 90 r 
' tan ~' ( y vs9 /|z vs 9 l) ! - 90 ° - 


YANG <90 


4. 4. 2. 2 Support parameters for the soft X-ray experiment (MA-048): The 

additional support parameters for the soft X-ray experiment are identical to 
those of the extreme-UV survey experiment. The only exception is that MA-048 
does not require any field-of-view computations. Refer to section 4. 4. 2. 4 for 
a detailed definition of the parameters . 

4. 4. 2. 3 Support parameters for the UV absorption experiment (MA-059): In 

the nominal mode of operation, only one parameter (THETA) will be computed in 
support of the UV absorption experiment. In addition, other parameters will be 
obtained from the reconstruction of the relative trajectory (Apollo vs. Soyuz) 
with use of the HOPE orbit determination program, on the basis of available 
relative tracking data. Special purpose utility programs will be used to merge 
THETA with the HOPE derived parameters d and Y , where d Is the relative distance 
between the CSM and the Soyuz and y is the angle between the relative position 
vector and the CSM velocity vector. Figure 4-4 illustrates these parameters. 


THETA 

The parameter THETA is computed as follows. 


-90° 1 THETA i 90° 


4. 4. 2. 4 Support parameters for the extreme-UV survey experiment (MA-083): 
The additional support parameters that are required to support the extreme-UV 
survey are specified below. These parameters also pertain to the MA-048 and 
HA-088 experiments. 

Since the MA-083 and MA -088 experiments require field-of-view computations, 
it is necessary to define the navigation base to nominal instrument base 
coordinate transformation F. The required transformation is defined as follows. 
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— T T T 
x = RDG 
x 


THETA = 90° - cos" 1 


x *v n 
x vl 


'vl' 



Instrument 

LOS 



Figure 4-4.- MA-059 support parameters. 


F = I y 


where 


z = (cos 6, sin 0 sin $ , -sin 0 cos <J>)‘ 
y = z x (l, 0, 0) T 


y “ y/|y 


x = y x z 


The nominal values for 0 and 4> are 90.0° and 37-75°, respectively. 


Geomagnetic position 

The geomagnetic position vector (designated Xll, Yll, and 7.1l) is computed 
as follows . 


r vll = B r vU 


The angle between the Earth-centered spacecraft position vector and the 
Earth-Run vector, SCRA, is computed in the following manner. 


0 = cos 


r sl‘ r vl 


l r sli ' r vll 


n=r. x r ,/ir . x r . | 
si vl 1 si vl 1 


s = n X r =1 , where r sl = ? sl /|r sl , 


0 if s*r . > 0 
vl — 


-8 if s-r , < 0 
vl 


• -l8o° < RCSA < l80° 


T**34«** • ' 


T 


ALF2 and DLT2 


The right ascension and declination of the vehicle in the ecliptic mean of 
1950*0 system is specified as follows. 


r „ = E r .. 
v2 vl 


ALF2 = tan" ! (y v2 /x v2 ) ; 0° <_ ALF2 < 360° 


DLT2 = sin-^z // r f ) ; -90° <_ DLT2 <_ 90° 


LOS and FOV orientation angles 


The computation of the instrument line-of-sight and field-of-view orientation 
angles is specified below for the following coordinate systems. 


1. Mean of 1950.0 

2. Ecliptic mean of 1950.0 

3. Up, east, north 


Note that the instrument alinement corrections are applied to the part 1 computa- 
tions and the FOV computations . 


1 . 


nl 


ulj 


L1RHA1 = tan _1 (y / x ' 
^nr nl' 


0 < L1FHA1 < 360 


L1DEC1 = sin _1 (z .) 


nl' 


-90° < LIDFjCI < 90° 


L1AE1 = L1HHA1 + AA1 
L1AD1 - L1DEC1 + ADI 


0° < L1AR1 < 360° 


-90° < L1AD1 < 90° 


where AA1 and ADI are the alinement corrections. 


For the field-of-view computations, replace the [0, 0, l] vector with the 

MA-083 field-of-view pointing vectors (x. 0 , x RQ , x nQ , x no ) or the corresponding 

Ay -tty uy uy 

MA-088 field-of-view pointing vectors , and compute 


F3P.HA1, F3DCA1 - right ascension and declination for the "A" 
quadrant for MA-083 
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F3RHB1, F3DCB1 - right ascension and declination for the ’’B” 
quadrant for MA-G83 


F8RHD1, F8DCD1 - right ascension and declination for the "D" 
quadra,:. ; for MA-088 

with the alinement corrections included in the computations. The essential 
difference "between the MA-O83 FOV and the MA-0R8 FOV is the angular width, 

2 degrees versus 15 degrees, respectively. 

0° <_ L1 P HA2 < 360° 

-90° <_ L1DEC2 <_ 90° 

0° £ LlA7i < 360° 

-90° <_ L1EL <_ 90° 

Vehicle velocity orientation angles 

The computation of the vehicle velocity orientation angles for the following 
coordinate systems are described below, 

1. Mean of 1950.0 

VRHA1 = tan -1 (y vl /x vl ) 

VDFC1 = sin’ 1 (7. ! ) 

vl 1 vl 1 

2. Ecliotic mean of 1950.0 


0° i VRHA2 < 360° 
-90° 1. VDEC2 <_ 90° 


v = E v 
v2 vl 


VRHA2 = tan’ 1 (v Jy_ ) 
v2 v2 


VDEC2 = sin’ 1 ( a _/l v !) 

v2 1 v2 1 ' 


0° £ VP.HA1 £ 360° 
-90° £ VDECi £ 90° 


x _ = E x .. 
n2 nl 


L1RHA2 = tan -1 (y Jx r 
n2 n^ 

L1DEC2 = sin _1 (z _) 
n2 


x n = CWT x , 
nlO nl 


L1AZ = tan _1 (y nl0 / 2nl0 ) 


L1EL = sin Kx 

nlO 


U6 




L0NS2 


The computation of the longitude of the Sun in the ecliptic mean of 1950.0 
coordinate system is as follows. 


r 0 = E r n 
s2 si 


L0NS2 = tan ‘(y^/x g) 


; 0° 1 LONS 2 < 360° 


YLOS 


The 
sight is 


angle between the vehicle velocity vector and the instrument line-of- 
computed as follows. 


YLOS « cos" 1 



; 0 < YLOS 1 180° 


ESLOS AND ESLOSS 


The angle between the Earth-Sun line and the instrument line-of-sight is 
computed in the following manner. 


ESLOS = cos" 1 




|r sll ‘ X nl ‘ 


ESLOSS = 180° -ESLOS 


0° < ESLOS < i8n° 


h.h.2.5 Support parameters for the helium glow experiment ( M A-Qfl8): The 

support parameters for the helium glow experiment are identical to those for 
the UV survey experiment (MA-083). .i detailed definition of the support 
parameters is presented in the preceding section. 



4. 4. 2. 6 Support parameters for the light flash experiment (MA-106): The 

additional parameters that are required to support the light flash experiment are 
presented below. 

ALPH10 , BETA10, and PHI10 

The orientation angles for the vehicle body axes with respect to the UEN 
coordinate system are computed as follows . 


_ m m m _ 

= GDR T ire x 10 


X 7 ~ YX 10 ^ y ij^ X 10 


ALPH10 = tan^y^/y^); °° 1 ALPH10 < 360° 

BETA10 = cos~ 1 (y 11 ); 0° £ BETA1G £ l8n° 

PHI10 = tan -1 (-y 21 /-y 31 ); 0° £ PHI10 < 360° 

where ALPIilO, BETA10, and PHI10 are illustrated in figure 4-5. 

BMAG and LMAG 

The magnitudes of the intensity of the Earth's magnetic field and the "L-shell" 
radius are computed as a function of vehicle position with the use of the NASA/ 

GSFC subroutines described in reference 9. The actual spherical harmonic 
coefficients that will be used to compute the geomagnetic field are designated as 
GSFC /65 in the above-referenced document. This model has been used in past 
programs to support experiments such as S-230 (Skylab). 

4. 4. 2. 7 Support parameters for the biostack experiment (MA-107): The 

additional parameters that are required to support the biostack experiment are 
the magnitude of the geomagnetic field intensity and the "L-shell" radius. The 
compuatations of these two parameters are discussed in the preceding section. 

4. 4. 2. 8 Support parameters for the Earth observations experiment (MA-136): 

The additional parameters required to support the Earth observations experiment 
are similar to those provided for Apollo except that the reference body is the 
Earth instead of the Moon. Also, the parameter list will be expanded to provide 
some parameters that were not computed for Apollo support. The required 
formulations for all parameters are presented below. 

The transformation F 1 for MA-136 is defined as 

F' = VAF 

[Fl 


F 



where 


7 , - (-cos B, -sin 0 sin «j> , sin 0 cos ) ‘ 

* m 

v = z x (l, 0, 0) ' 

y = y/|y| 

x = y x z 


and where 6 = 57-5° and <j> = 58.9°. 

The transformation A r orients the camera x-axis to within +4 5° of the vehicle 
velocity vector; V is specified "below. 


V = 


cos(-niT/2) sin(-nir/2) 0 
-sin(-nir/2) cos(-nir/2) 0 
0 0 1 


where 


n = k if — tt/4 + kiT/2 < 0 < it/ 4 + kn/2 for k = 0, 1, 2, and 3 


i 


and where 

0 = tan _;L (v y0 /v x0 ) 
\q = AFGDR ~ vl 


See figure 4-6 for ar illustration of the nominal instrument base coordinate system. 
Principal point and FOV footprint 

The computation of the principal point and the field-of-view intersection 
points is presented below. 

T T T T- 

- WT R D G F' x i9 ; i = A, B, C, D, n 

where A, B, C, and D designate the F^V and n designates the nrincipal point. Goromte 
W ? a4’ ? b 4’ ? C4’ 811(1 ? D4 algorithm 1. 
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Figure 4-6 MA-136 camera coordinate system (CM right side window). 
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MT\ = tan -1 1 z^/px ^ 2 + y^ 2 )^ 2 (l - f) 2 J J ; - 90 ° < LAT ± <_ 9 0° 
LON. = tarTKy^/x.^); 0° < LON. < 360°' 

Algorithm X 
Let 

r = principal point or FOV intersection vector with respect to origin 
? v = position of the vehicle 


x. = instrument pointing vector (unit vector) 
d = magnitude of instrument pointing vector 

where all vectors are with respect to the geographic inertial coordinate system. 
Nov 


r = r + dx. 
v 1 

2 2 2 

f(?) 

a a b 

Given an initial estimate d, denoted by d = r cos a - */ r 2 - r 2 sin 2 « 

o v ™ v 

where 


r = geodetic radius of the Earth at ths subvehicle point 
r^ = magnitude of r^ 

a = angle between x. and -r 
1 v 


compute 


r = r + d x. 
o v 01 


f (r ) = k 
o o 
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To estimate d 


r\ rr» A _ 

f(r) - f( r ) = 1 - k = -Pr- <5d = m6d 
o o ou 

o 

d = d + <5d 
o 

where 



3r 

o — 

3d = X i 

5d = (l - l'- 0 )/m 

Repeat the process until |l - k^j < e. Set r = r^ . 

Slant range (SR) 

The slant range SR is a by-product of algorithm 1; that is, 

SR = d 

Scale factor (SF ) 

The scale factor is computed as follows. 

SF = f/h 

where h is the geodetic altitude of the vehicle and f is the focal length of 
the camera lens . 

Focal length (FL) 

The camera focal length is available by card input as a function of time. 

Altitude rate (ALTR) and horizontal velocity (HV) 

The altitude rate and the horizontal velocity with respect to the UEN 
coordinate system associated with the principal point are defined as follows. 
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v vlO 
HV = 

ALTP 


= CWT v . 

P vl 

(W + v zl0 2 ) 

= v xlO 


l/s 


where C designates the UEN system associated with the principal point. Figure 4-7 
illustrates these parameters. 

TILT and TILTAZ 

These parameters are computed as follows. 


_ T T T T 

x nlO = <VWTrVg T F ,T 


0 

0 

-1 


TILT = cos"' 1 (-x nl0 ) ; 0° £ TILT < 180° 


TILTAZ = tan -1 (y ../z 0° < TILTAZ < 360< 

nlO nlO 5 — 


where designates the UKN system associated with the vehicle, ™igure 4-8 
illustrates these parameters. 

uELF_ and f'AZP 

The Sun azimuth and Sun elevation with resrect to the principal point are 
specified as follows. 


r , = c W r . 

slQ p si 


"pslh r sl^ r plf 

. -1 


SELP = sin -1 (x '_/ r . n ) ; -90° < SELP < 90 c 
psiu psiu — — 


1A7.P = n 0 < PAKP < 360° 


where C designates the UEN system associated with the principal point. See 
P 

figure 4-9, which illustrates the output parameters. 
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Up(X) 



V - inertial velocity of the vehicle 


Figure 4-7.- ALTR and HV parameters. 



\ Optical axis 
* (-Z camera axis ) 


SVP - subvehicle point 


Figure 4-8 Tilt and tilt azimuth. 
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LATS and LONS 


The geodetic latitude and longitude of the subsolar point are computed as 
follows . 


r s4 = WT ? sl 

LATS = tan" 1 | )* 2 (l “ f) 2 J '•> 

LONS = tan -1 (y slt /x sl4 ) j 0° £ LONS < 360° 


-90° < LATS < 90° 


These parameters are illustrated in figure 4-10 . 
EMISS, NDA, and SWING angles 

These parameters are specified as follows. 



rn m 

t? trjppw J ]-.rr 


-1 

0 

0 


EMISS = cos -1 (-z /|r u9 |) ; 0° £ EMISS £ l80° 


— T T 

r n«rtr*PT rC T 

n° i 


0 

0 

1 


IDA = tan" 1 (v Jx .)*, 0° < IDA < 36o° 
n u n« 

r i 

-i 

^ rn n rn 

r = T’T.pp T \r0 ' 0 

u9 v 

' 0 

FL'INO = tan” 1 (x _/v ) : Q° < SWP!0 < 3fn c 

u n * u° — 


where and designate the UEN systems associated with the principal point 

and the vehicle, respectively. Figures 4-11 through 4-13 depict these parameters. 




Optical axis 
(-Z camera) 



Y - camera 


PP ** principal point 


Figure 4**11 Emission angle. 
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Optical axis 
(-Z camera) 

-Up (SVP) 

^ X - camera 

Y - camera 



SVP - subvehicle point 


Figure 4-13 Swing angle. 
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PHASE 


The Sun phase angle is defined below. 



WT r 


si 


r ps4 *p4 



0° < PHASE < l80° 


Figure 4-l4 illustrates this parameter. 

XTILT, YTILT , and HEAD angles 

The specification of the XTILT, YTILT and HEAD angles is presented below. 
Refer to figure 4-15, which depicts these parameters. 

— TIT T 

X xlO = °v WTR D G ' F 

YTILT = sin -1 (x Q ); -90° <_ YTILT <_ 90° 

HEAD = tan -1 1C> ) 1 °° 1 HEAD < 360° 

_ T> T 1 T 1 T 

X yl0 = C V WTR‘D G F» 

XTILT = sin - 90 ° £ t jj jT hq 0 

where C_^ designates the UEN system associated with the vehicle. 

Ham era orientation angles 

The camera orientation (PHI, KAPPA, and OMEGA) angles with respect to the 
UEW coordinate system are comnuted as follows • 


0 

1 

0 


1 

0 

0 
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Up 00 



SVP - subvehicle point 

Figure 4-15.- XTILT, YTILT, and HEAD angles. 
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T T T 

M = F'GDRT WC 

v 


0 0 1 
10 0 
0 10 



PHI = tan -1 (m 13 /m 33 ) ; 0° <_ PHI < 360° 
KAPPA = sin" 1 (-m 23 ) ; -90° <_ KAlPPA _< 90° 
OMEGA = tan _1 (m 21 /m 22 ) ; 0° <_ OMEGA < 360° 


where C designates the UEN system associated with the vehicle, 
v 

L08X, 10SY , and LOSZ 

The direction cosines for the Vector from the spacecraft to the principal 
intersection point in the geographic system are computed in the following manner. 


L0SX 
L0SY 
LOSZ 

where x, y, and z are unit vectors in the geographic coordinate system. 
GEOCAM and GEOUEN 

The geographic to camera system transformation and the UEN (vehicle) to 
camera axis transformation are specified as follows. 

GC = E'ODRtV 

T T T 

LH = F'GDRT ¥ C 


ALPHA angle 

The A.LPHA angle, depicted by figure U-lb, is computed in the following 
manner . 


- x n H‘ x 


= x i, v 


nh 


- x , z 
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r psl* = W r sl - r p 4 
X n4 = 

r ps4 ~ r ps4 // | r ps i(l 
n ps4 " ( r ps 4 * x n4^ 
n ps4 “ n p S 4 / l n p S 4l 


si, = C 


t = s. 


ALPHA = cos 


( n ps4‘ s 4 ) n P s4 

cos 1 [(-X n4 .tJ/jt| ] ; 0° £ ALPHA < 180° 


Sur face arc length (AL) 

surface distance between nadiy nr ia +-U 

following manner. &nd the P^ncipal point is computed i, 


r v4 ~ ^v4 ^ I r \ 


r p4 r p 4 / ' r p4' 
0 = 

AL = r( A , tj>) 9 


''here r(j, ij>) is the geodetic radius at thf* n a- 

length computation is illustrated in £g^e “w! T ’° lnt ' The surface 

Forward overlap ratio ( DVR ) 


The forward 


overlap ratio is computed as follows. 
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Earth 

center 



Figure 4-17 Surface arc length . 
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•V'V * ' r d4 (t i ,/ l ? aU (t i ) i 

n(ti> = r ck (t.) x r al( (t.) 
n(t i ) = n(t i )/|n(t i ) | 

5 = ? d. ( Vl ) -[ :; c4 (t i+3 ),n(t i'l ntt i ) 

Compute the arc length between s and r^(t^), AL s » and the arc length between 
r , (t.) and 7,i (t.), AL Compute the overlap ratio as follows. 

C4 X 0.4 1 cd. 

OVR = AL /AL. 

s cd 

See figure H-l8. 

DELTAT 

The time interval between adjacent photographs is 
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5.0 EXPERIMENT SUPPORT PROGRAMS 


5.1 Design Philosophy 

There are several major programs and a large number of utility programs that 
will be involved in the preprocessing and processing of ASTP mission data in order 
to produce requested experiment support products. This section lists each of the 
major programs and provides a brief description of the functions and capabilities 
of each. Further information concerning any of these programs can be obtained 
through the referenced documents or by contacting the Mathematical Physics Branch. 

In formulating the ASTP experiment support system, the major design philosophy 
was to maximize the use of existing software. As a result, only one truly new 
program, the Apollo-Soyuz Experiment Parameter (ASEP) program, was designed; and 
even that program made extensive use of formulation from its lunar orbit predecessor, 
the Apollo Photo Evaluation (APE) program. All other programs used in the processing 
represent existing software. Only minor modifications to existing programs were 
required in order to accommodate ASTP data formats. In adhering to this philosophy, 
several elements of the ASTP system were designed so as to minimize the programing 
impact rather than maximize the efficiency and utility of the resulting design. 

Figure 5-1 shows the experiment support system design from a program standpoint and, 
as such, is a simplification of figure 3-1. In designing future systems, it is hoped 
that we may be able to more closely follow the design suggested in figure 5-2. The 
basic functions of this proposed system are identical to those of the ASTP system, 
but many of the functions have been combined and/or modified to facilitate the 
handling of large volumes of data. The differences between the proposed system and 
the ASTP system are enumerated below. 

1. The T/M preprocessing functions, performed in two steps by the GPREP and 
GEDIT programs for ASTP, would be combined into a single processor. The use of 
an online processor similar to the RTCC system, with DTV and user interface capa- 
bilities, could reduce the t/m data preprocessing operation (refers only to that 
portion of the T/M data reduction directly related to ephemeris, pointing, and 
field-of-view data) for a mission such as ASTP from an effort requiring several 
months of offline processing to one that could be completed in a day or less once 
all raw data become available at JSC. 

2. The trajectory generation function would provide only state vectors and the 
associated physical constants required for their propagation by other integrators. 

In the ASTP system, the trajectory tape output by the trajectory generation program 
(HOPE) eventually becomes the driver for the program that generates the experiment 
support data (ASEP). Thus, all desired output times must be known prior to 
generating the trajectory tape. As a result, the raw T/M preprocessing and 
trajectory generation must proceed serially. In the proposed system, the trajectory 
integration package would be contained within the experiment support program and 
the trajectory generation would be accomplished in the final step on the basis of 
the state vectors and constants provided. Thus , the raw T/M preprocessing and 
trajectory determination functions could proceed independently and in parallel. 

3. The experiment support programs, as mentioned above, would be modified to 
include a trajectory integration package. 
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Figure 5-1.- AS TP experiment support system software. 
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Figure 5-2.- Optimized experiment support system design. 
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It is believed that utilization of a system with these design characteristics 
could reduce production time for a mission such as ASTP from the current estimate 
of 4 to 5 months to a maximum of 1 month. Manpower requirements could be 
correspondingly reduced by 60 to 80 percent. 


5.2 Preprocessing Programs 

5.2.1 GPREP . - GPREP is ah attitude data preprocessing program that strips 
time and CSM gimbal angles from a tape containing data from the CMC telemetry 
downlist and reformats these data in a form compatible with the experiment support 
processing. 

5.2.2 GEDIT . - The GEDIT program operates on tapes output by GPREP. The 
purpose of the program is to provide a clean attitude data tape to be used in the 
experiment support data production. The raw attitude data are edited by eliminating 
duplicate points, discarding all bad data, and, where possible, filling data gaps 
with the use of a linear interpolation routine. (The use of a higher order routine 
was investigated, but the linear routine proved to be best suited for the job 
because of the linear nature of the orientation angles.) The program is also 

used to merge various segment;,, of edited data to form master attitude data tapes 
and finally to merge these with HOPE-generated ephemeris tapes to form the merged 
trajectory and gimbal tapes that are used as direct input to the experiment 
support (ASEP) program. 

5.2.3 LOVE .- The LOVE program is used to preprocess spacecraft tracking data. 
It has the capability of generating master tracking data tapes from raw data in 
any of the formats currently provided by the STDN network, GSFC, and JPL. 


5.3 Processing Programs 

5.3.1 HOPE . ~ The Houston Operations Predictor/Estimator (HOPE) - references 3, 

4 , and 5 - is a two-vehicle double-precision orbit determination/prediction program. 

It has the capability of solving for the inertial and/or relative orbits of one or 
two orbiting vehicles with the use of a variety of observation types, including 
ground-based, vehicle-vehicle , and vehicle-ground observations. The orbit determination 
is accomplished by means of a weighted-least-squares differential correction 
technique. The program also has stand-alone integration capability for use as a 
trajectory prediction tool, as well as dummy data generation capability. The user 

has ready access to all program constants by means of simplified card and/or tape 
input . 

5.3.2 ASEP .- The Apollo-Soyuz Experiment Parameter (ASEP) program acts as an 
executive routine for the generation of spacecraft ephemeris and orientation data 
from the merged trajectory and gimbal tapes (see sec. 5 . 2.2 and fig. 3-1), in 
addition to calculating all special experiment-related parameters requested by 
various Pi's. The program outputs all PI experiment support data (microfilm, 
tabulator listings, and experiment support data tapes). Details of the ASEP design 
can be found in section 6 . 



6.0 EXPERIMENT SUPPORT FROGRAM DESIGN 


The purpose of this section is to specify the design of the Apollo-Soyuz 
Experiment Parameter (ASEP) program that will provide ephemeris, attitude, 
pointing, and field-of-view parameters as defined in section 4.0. An overview 
of the program is presented in section 6.1. The details of the program design 
are discussed in section 6.2, whereas the required subroutines are described in 
section 6.3. 


6.1 Program Structure 

To provide a better understanding of the ASEP program, a brief discussion of 
the major program functions is presented in this section. The major functions of 
the program can be divided into three categories: input processing, support 

parameter computations, and output processing. Figure 6-1 provides an overview of 
these functions. 

The first step is to read and process the input control cards and the mission 
data cards. Next, the input data tape is searched until the processing start time 
(input by card) is found. Now the computation of the support parameters begins. 

On the basis of the trajectory and attitude data that are read from tape, the 
nominal support parameters are computed as specified in section 4.4.1. Next, the 
experiment on-off times for MA-007 are tested to determine if the parameters 
peculiar to that experiment are required. If they are, the parameters are 
computed as specified in section 4. 4. 2.1. If not, then the on-off times for the 
next experiment are tested. This process continues until all required experiment 
specific parameters are computed. 

At this point, all the output parameters are formatted and output on tape, 
microfilm, and listing as required. When this is completed, the next data 
record is read from the input data tape and the process is repeated until all 
the data records are processed or the processing stop time (input by card) is 
reached. The final operation is to write an end-of-file on the output tape. 

6.2 Detailed Design 

This section specifies in detail the various processing elements of the ASEP 
program: input - tape and card, computations (nominal and experiment specific) - 

and output - tape, microfilm, and listing. 

6.2.1 Input . - 

6. 2.1.1 Tape: The input tape that is produced by the GEDIT program (see 

sec. 3.2,l) provides time-ordered ephemeris, attitude, and transformation data 
from which all the experiment support parameters can be computed. Table 6-1 
specifies the format of this tape. 

6. 2. 4. 2 Card: Card inputs are required for the specification of ASEP program 

options and for loading ASTP mission and experiment specific data. Table 6-II 
specifies the name, type, dimension, and a description of the function of each of 
the input variables. The Univac 1108 Exec II NAMELIST option will be used to 

load the data. 
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Figure 6-1 ASEP processing overview. 





















TABLE 6-1.- ASEP INPUT DATA TAPE FORMAT 




(a) Initial 

Word 

15 

Type 

1 

BTIME 

DP 

2 

BTIME 

DP 

3 

BTIME 

DP 

4 

BTIME 

DP 

5 

BTIME 

DP 

6 

BTIME 

DP 

7 

GET 

DP 

8 

GHAE 

DP 

9 

GHAL 

DP 

10-20 

TITLE 

HOL 

21-23 

N/A 

SP 


(b) 

Data record (units: 

Word 

ID 

Type 

1 

TAET 

DP 

2-4 

PECL 

DP 

5-7 

VECL 

DP 

8-10 

AECL 

DP 

11-13 

PTOD 

DP 

14-16 

VTOD 

DP 

17-19 

ATOD 

DP 

20-22 

PGE0 

DP 

23-25 

VGEO 

DP 

26-28 

AGEO 

DP 


record 

Description 
Year (mod. 1900) 

Month 

Day 

> Base date 

Hour 

Minute 

Second ^ 

Minutes from launch of base time 

Right ascension of Greenwich midnight day 
of epoch 

Right ascension of Greenwich midnight day 
of launch 

66 alphanumeric characters from HOPE title 
card 

Not used 

km, sec, and deg) 

Description 
Time (min) from BTIME 

Vehicle position - ecliptic mean of 1950.0 
Vehicle velocity - ecliptic mean of 1950.0 
Vehicle acceleration - ecliptic mean of 1950.0 
Vehicle position - true of date 
Vehicle velocity - true of date 
Vehicle acceleration - true of date 
Vehicle position - geographic inertial 
Vehicle velocity - geographic inertial 
Vehicle acceleration - geographic inertial 
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TABLE 

6-1.- ASEP 

INPUT DATA 

TAPE FORMAT - Concluded 


(b) Data 

record (units: km, 

sec, and deg) - Concluded 

Word 

ID 

Type 


Description 

29-31 

PM50 

DP 


Vehicle position - mean of 1950.0 

32-34 

VM50 

DP 


Vehicle velocity - mean of 1950.0 

35-37 

AM 50 

DP 


Vehicle acceleration - mean of 1950.0 

38-40 

PROT 

DP 


Vehicle position - geographic rotating 

41-43 

VROT 

DP 


Vehicle velocity - geographic rotating 

44-46 

AROT 

DP 


Vehicle acceleration - geographic rotating 

47-55 

EMAT 

DP 


Mean of 1950.0 to ecliptic mean of 1950.0 
matrix 

56-64 

TMAT 

DP 


Mean of 1950.0 to true of date matrix 

65-73 

WTMAT 

DP 


Mean of 1950.0 to geographic inertial matrix 

74-76 

SUNP 

DP 


Sun position - mean of 1950.0 

77-79 

SUNV 

DP 


Sun velocity - mean of 1950.0 

80-82 

CDUX, 

CDUY, 
CDUZ 

DP 


Gimbal angles 

83-85 

MUNP 

DP 


Moon position - mean of 1950.0 

86-88 

MUNV 

DP 


Moon velocity - mean of 1950.0 

89 

GIMF 

DP 


Gimbal status flag 

90-91 

N/A 

DP 


Not used 



(c) 

Terminating record 

Word 


Type 


Value 

1-91 


DP 


0.00 

92 


SP 

T9 

EOFAAA (where A = blank) 



TABLE 6-II.- INPUT CARDS FOR ASEP 


Variable 


name 

Type 

Dimension 

Function 

DRIFTX 

DP 

(4,4) 

IMU X-drift (deg/hr) as a function of time 
(hr, min, sec) 

DRIFTY 

DP 

(4,4) 

IMU Y-drift (deg/hr) as a function of time 
(hr, min, sec) 

DRIFTZ 

DP 

(4,4) 

IMU Z-drift (deg/hr) as a function of time 
(hr, min, sec) 

EXP007 

DP 

(3,2,10) 

Operation periods for MA-007; start time 
(hr, min, sec) and stop time (hr, min, sec) 

EXP048 

DP 

(3,2,10) 

Operation periods for MA-048 ; start time 
(hr, min, sec) and stop time (hr, min, sec) 

EXP059 

DP 

(3,2,10) 

Operation periods for MA-059; start time 
(hr, min, sec) and stop time (hr, min, sec) 

EXP083 

DP 

(3,2,10) 

Operation periods for MA-083; start time 
(hr, min, sec) and stop time (hr, min, sec) 

EXP088 

DP 

(3,2,10) 

Operation periods for MA-088; start time 
(hr, min, sec) and stop time (hr, min, sec) 

EXP106 

DP 

(3,2,10) 

Operation periods for MA-106; start time 
(hr, min, sec) and stop time (hr, min, sec) 

EXP136 

DP 

(3,2,10) 

Operation periods for MA-136; start time 
(hr, min, sec) and stop time (hr, min, sec) 

FL 

DP 

(4,11) 

Lens focal length (mm) as a function of tim^ 
(hr, min, sec) 

GTBIAS 

DP 

(3,2,6) 

g.e.t. time bias (hr, min, sec) as a function 
of time (hr, min, sec) 

LAUNCH 

DP 

(6) 

Launch time (yr (1900), month, day, hr, min, sec 

L1ADC1 

DP 

(4,6) 

Alinement correction (deg) for declination for 
MA-083 as a function of time (hr, min, sec) 

L1ARA1 

DP 

(4,6) 

Alinement correction (deg) for right ascension 
for MA-083 as a function of time (hr, min, sec) 

MODPRT 

DP 

(1) 

Control flag for output print frequency 

PET 

DP 

(3,2,5) 

PET time bias (hr, min, sec) as a function of 
time (hr, min, sec) 

REV 

DP 

(1) 

Rev., number referenced to TSTART 


r 
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TABLE 6-II INPUT CARDS FOR ASEP - Concluded 


Variable 


name 

Type 

Dimension 

Function 

RF 

DP 

(3,4,11) 

REFSMMAT matrix by row as a function of time 
(hr, min, sec) 

TSTART 

DP 

(3) 

Start time of run (hr, min, sec) 

TSTOP 

DP 

(3) 

Stop time of run (hr, min, sec) 
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Card input (NAMELIST) will also be used to override nominal program and experi- 
ment constants. Tables 6-III and 6- IV define -all required program constants, the 
variable names, and the nominal values of the constants. 

6.2.2 Nominal support parameters .- For each time point within the interval 
specified by TSTART and TSTOP, the nominal support parameters will be computed. 

Figure 6-2 illustrates the straightforward procedure by which all the required 
parameters are computed. Note that besides the nominal computations, matrices C, 

G, D, and R are precomputed for later use in the computation of the experiment- 
specific parameters. 

6.2.3 Experiment-specific support parameters .- After the nominal parameters 
are computed, the experiment on-off times are tested to determine if the experiment- 
specific parameters need to be calculated. The sequential testing procedure and 
associated computations are illustrated in figure 6-3. Note that if any one 

of the three experiments MA-0U8, MA-083, and MA-088 are operating, then the parameter 
set, common to all three experiments, is computed. Parameter statistics are also 
computed at this point. 

6.2 . k Output . - The forms of the output data include tape, microfilm, and 
listing. These outputs are discussed in detail in section 7 and in appendix B. 

6.3 Subroutine Descriptions 

This section describes the function of the various subroutines used by the 
ASEP program. These routines are described in tables 6-V and 6-VI. 


TABLE 6- I I I.- ASEP PROGRAM CONSTANTS 


Parameter 

Nominal Value 

Description 

AE 

6.378166000D3 

Semimajor axis of the Earth's refer 
ence ellipsoid (km) 

BE 

6.356784287D3 

Semiminor axis of the Earth's 
reference ellipsoid (km) 

CRAD 

5 . 72957795130823209D1 

Radians to degrees 

EMU 

3. 986032D5 

3 2 

Mu of the Earth (km / sec ) 

ER2KM 

6 . 378165D3 

Earth radii to kilometers 

OMEGE 

4.3752690880008D-3 

Angular velocity of the Earth 
( rad/min) 

RADS 

6 . 96495618D5 

Radius of the Sun (km) 

REM 

6.367475142D3 

Mean Earth radius (km) 

RS 

1.49597900D8 

AU (km) 

TWOPI 

6. 28318530717958648D0 

2tt 
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TABLE 6- IV.- ASEP EXPERIMENT CONSTANTS 


Parameter 

Nominal value 

Description 

CLTMAG 

11.5 

Colatitude of the Earth's magnetic 
North Pole (deg) 

F0V083 

1.0 

Field-of-view for MA-083 (deg) 

F0V088 

7.5 

Field-of-view for MA-088 (deg) 

FOVL36 

53.0 

Field-of-view for MA-136 (mm) 

LONMAG 

-69.5 

Longitude of the Earth's magnetic North 
Pole (deg) 

MA007 

*3x3 

MA-007 misalinement matrix 

MA083 

*3x3 

MA-083 misalinement matrix 

MAI 36 

*3x3 

MA-136 misalinement matrix 

PHI007 

58.9 

^ body orientation angle for MA-007 (deg) 

PHI083 

37.75 

<f>, body orientation angle for MA-083 (deg) 

PHI 136 

58.9 

<J>, body orientation angle for MA-136 (deg) 

THT007 

57.5 

6, body orientation angle for MA-007 (deg) 

THT083 

90.0 

0, body orientation angle for MA-083 (deg) 

THT136 

57.5 

0, body orientation angle for MA-136 (deg) 


84 



TABLE 6-V.“ UTILITY SUBROUTINES 


Subroutine 


name 

Function 

DPt-M 

Multiplies matrices 

MATRN 

Transposes matrices 

ROTX 

Computes single-axis rotation 
matrix 

UNIT 

Vector operations - unitization, 
magnitude, dot product, and cross 
product 




TABLE 6- VI.- COMPUTATIONAL SUBROUTINES 


Subroutine 

name 

ADBARV 

CARMEL 

DRIFT 

FOV 

GI2GE0 

GI2UEN 

IM2NAV 

INTEG 

INVAR 

KEPLER 


Function 

Computes spherical polar elements 
(see sec. 4. 4. 1.2) 

One of six NASA/GSFC subroutines 
used to compute B&L coordinates 

Computes the D matrix (see 
sec. 4.3.2) 

Computes the MA-136 LOS and FOV 
intersection points 

Computes the B matrix (see sec. 4.3.10) 

Computes the C matrix (see sec. 4.3.9) 

Computes the G matrix (see sec. 4.3.3) 

One of six NASA/GSFC subroutines 
used to compute B&L coordinates 

One of six NASA/GSFC subroutines 
used to compute B&L coordinates 

Computes Keplerian elements (see 
sec. 4. 4. 1.2) 



TABLE 6-VI COMPUTATIONAL SUBROUTINES - Concluded 


Subroutine 

name 

LINES 

MA007 

MAO 48 

MAO 59 

MA106 

MA136 

NEWMAG 

START 

SUBVEH 

SUNANG 


Function 

One of six NASA/GSFC subroutines 
used to compute B&L coordinates 

Computes the support parameters 
for the MA-007 experiment (see 
sec. 4.4.2. 1) 

Computes the support parameters 
for the MA-048/MA-083/MA-088 
experiments (see sec. 4. 4. 2. 4) 

Computes the support parameter 
for the MA-05S experiment (see 
sec. 4.4. 2. 3) 

Computes the support parameters 
for the MA-106 experiment (see 
sec. 4.4. 2. 6) 

Computes the support parameters 
for the MA-136 experiment (see 
sec. 4. 4. 2. 8) 

One of six NASA/GSFC subroutines 
used to compute B&L coordinates 

One of six NASA/GSFC subroutines 
used to compute B&L coordinates 

Computes the geodetic altitude, 
latitude, and longitude of the 
vehicle (see sec. 4.4. 1.4) 

Computes the Sun azimuth and ele- 
vation angles (see sec. 4. 4. 1.4) 
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7.0 EXPERIMENT SUPPORT OUTPUT DATA FORMATS 


7 . 1 Output Tape 

The master support data tape consists of 7 records of input, constants, and 
format information, followed by the 2^1-word logical data records. Each data tape 
will be terminated with an end-of-file (EOF). 

7.1.1 Header records .- The seven information records that precede the data 
records include the following items. 

1. Input data - four records. These records include all namelist input data. 

2. Program constants - one record. This record includes all the physical 
constants and experiment-specific data that will be used to compute the support 
data. 


3. Data format - two records. These records specify the name and location 
on tape of all the support parameters. 

It should be noted that the seven header records appea~ -nly on the master data 
tapes retained at JSC. They are not a part of PI CC1 

7.1.2 Data records . - 2kl single-precision words. The format for the de/.a 
records is specified in table 7-1. The table lists the parameter name, the type 
of parameter, the parameter location, and the parameter definition. For cross- 
reference purposes , table 7-II lists the output parameters in alphabetical order, 
with the associated word .location. Note that the parameter units are kilometers, 
seconds, and degrees unless otherwise specified. 

7.2 Listing and Microfilm 

Figure 7-1 illustrates the output format for both listing and microfilm. All 
the parameter names correspond to those in table 7-1 « The frequency of the micro- 
film output will be the same as the tape output frequency, whereas the listing 
frequency will be some integer multiple of the tape output frequency as specified 
by card input (MODPRT) . 

7*3 Output Parameter Statistics 

As a quality assurance procedure, output parameter statistics will be computed 
for each support parameter. The statistics to be computed are as follows. 

Mean 

Second moment about the mean 

Third moment about the mean 

Fourth moment about the mean 
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TABLE 7-1.- EXPERIMENT SUPPORT DATA TAPE FORMAT 


Parameter name(s) 

Type of 
parameter 

Parameter 

location 

Parameter definition 

GMT 

Nominal 

1-6 

Greenwich mean time (yr (mod. 1900), month, 
day, hr, min, sec) 

OPFLAG 

Nominal 

7 

Code word to designate which experiments 
are operating 

AETH, AETM, AETS 

Nominal 

8-10 

Ground elapsed time (hr, min, sec) 

GETH, GETM, GETS 

Nominal 

11-13 

ASTP mission elapsed time (hr, rnin, sec) 

CTEH, CTEM, CTES 

Nominal 

14-16 

Onboard clock elapsed time (hr, min, sec) 

PETH, PETM, PETS 

Nominal 

17-19 

Phase elapsed time (hr, min, sec) 

XI, YI, Zl, XD1, YD1, ZD1 

Nominal 

20-25 

State vector - Cartesian elements in the mean 
of 1950.0 coordinate system 

ALF1, DLT1, BTA1, AZl, Rl, 
VI 

Nominal 

26-31 

State vector - ADBARV elements in the mean 
of 1950.0 coordinate system 

SMA1 , ECC1 , INC1, NODI, 
OMGl , TA1 

Nominal 

32-37 

State vector - the Keplerian elements in the 
mean of 1950.0 coordinate system 

X12, Y12, Z12 , XD12, 
YD12 , ZD12 

Nominal 

38-43 

State vector - Cartesian elements in the 
geographic rotating coordinate system 

ALF12, DLT12, BTA12, 
AZ12, Rl 2 , VI 2 

Nominal 

44-49 

State vector — ADBARV elements in the 
geographic rotating coordinate system 

SMA12, ECC12, ING12, 
NOD12, OMGl 2, TA12 

Nominal 

50-55 

State vector - Keplerian elements in the 
geographic rotating coordinate system 

X3 , Y3, Z3, XD3, YD3, ZD3 

Norn in al 

56-61 

State vector - Cartesian elements in the true 
of date coordinate system 



\o 

rv> 




Ell, E12, E13 


Nominal 


Parameter 

location Parameter definition 

62-67 State vector - ADBARV elements in the true of 

date coordinate system 

68-73 State vector - Keplerian elements in the true 

of date coordinate system 

74-79 State vector - Cartesian elements in the 

geographic inertial coordinate system 

80-85 State vector - ADBARV elements in the 

geographic inertial coordinate system 

86-91 State vector - Keplerian elements in the geographc 

coordinate system 

92-97 Solar Cartesian elements in the mean of 1950.0 

coordinate system 

98-100 Row 1 of REFSMMAT 

101-103 Gimbal angles 

104-106 Row 2 of REFSMMAT 

107-109 Geodetic altitude, latitude, and longitude 

110-112 Row 3 of REFSMMAT 

113-114 Solar azimuth and elevation wrt the subvehicle 

point 

115 Greenwich hour angle 

116-118 Row 1 of the T matrix (mean of 1950.0 to true 

of date) 

119-121 Row 1 of the E matrix (mean of 1950.0 to ecliptic 

mean of 1950.0) 




TABLE 7-1.- EXPERIMENT SUPPORT DATA TAPE FORMAT - Continued 
Type of Parameter 


Parameter name(s) 

parameter 

location 

T21 , T22, T23 

Nominal 

122-124 

E21, E22, E23 

Nominal 

125-127 

T31 , T32, T33 

Nominal 

128-130 

E31 , E32 , E33 

Nominal 

131-133 

GIMB 

Nominal 

134 

REV 

Nominal 

135 

SUNF 

MA-007 

136 

SUNH , SUNM, SUNS 

' MA-007 

137-139 

XANG, YANG 

MA-007 

140-141 

XI 1, Yll, Zll 

MA-048 , 
MA-083 , 
MA-088 

142-144 

SCSA 

MA-048 , 
MA-083, 
MA-088 

145 

ALF2, DLT2 

MA-048 , 
MA-083, 
MA-088 

146-147 

L1RHA1, L1DEC1 

MA-048, 

MA-083, 

MA-088 

148-149 

L1RH2, L1DEC2 

MA-048, 

MA-083, 

MA-088 

150-151 


■Parameter definition 
Row 2 of the T matrix 
Row 2 of the E matrix 
Row 3 of the T matrix 
Row 3 of the E matrix 
Gimbal status flag 
Rev. number 
Sunrise/sunset flag 
Sunrise/sunset time (hr, min, sec) 
x- angle and y-angle deviations 
Geomagnetic position of the vehicle 

Vehicle - Sun angle 


Right ascension and declination of the vehicle in 

the ecliptic mean of 1950.0 coordinate 

system 

Right ascension and declination of the LOS in the 
mean of 1950.0 coordinate system 

Right ascension and declination of the LOS in the 
ecliptic mean of 1950.0 coordinate system 



TABLE 7-1.- EXPERIMENT SUPPORT DATA TAPE FORMAT - Continued 




Type of 

Parameter 



Parameter 

name ( s ) 

parameter 

location 

Parameter definition 


L1AZ, L1EL 

MA-048 , 

152-153 

Azimuth and elevation of the LOS 




MA-083 , 






MA-088 




F3RHA1 , 

F3DCA1 

MA-048 , 

154-155 

Right ascension and declination of FOV A 

in 



MA-083, 


the mean of 1950.0 coordinate system for 




MA-088 


MA-083 


F3RHB1, 

F3DCB1 

MA-048, 

156-157 

Right ascension and declination of FOV B 

in 



MA-083 , 


the mean of 1950.0 coordinate system for 




MA-088 


MA-083 


F3RHC1, 

F3DCC1 

MA-048, 

158-159 

Right ascension and declination of FOV C 

in 



MA-083, 


the mean of 1950.0 coordinate system for 




MA-088 


MA-083 


F3RHD1 , 

F3DCD1 

MA-048 , 

160-161 

Right ascension and declination of FOV D 

in 



MA-083 , 


the mean of 1950.0 coordinate system 




MA-088 


for MA-083 


F8RHAI , 

F8DCA1 

MA-048 , 

162-163 

Right ascension and declination of FOV A 

in 



MA-083, 


the mean 1950.0 coordinate system for 




MA-088 


MA-088 


F8RHB1, 

F8DCB1 

MA-048, 

164-165 

Right ascension and declination of FOV B 

in 



MA-083 , 


the mean 1950.0 coordinate system for 




MA-088 


MA-088 


F8RHC1, 

F8DCC1 

MA-048, 

166-167 

Right ascension and declination of FOV C 

in 



MA-083, 


the mean 1950.0 coordinate system for 




MA-088 


MA-088 


F8RHDI, 

F8DCD1 

MA-048 , 

168-169 

Right ascension and declination of FOV D 

in 



MA-083 , 


the mean 1950.0 coordinate system for 




MA-088 


MA-088 




TABLE 7-1.- EXPERIMENT SUPPORT DATA TAPE FORMAT - Continued 


Parameter name(s) 

Type of 
parameter 

Parameter 

location 

Parameter definition 

L1ARA1 , L1ADC1 

MA-048 , 
MA-083, 
MA-088 

170-171 

Corrected right ascension and declination of the 

LOS in the mean of 1950.0 coordinate 

system 

VRHAi, VDEC1 

MA-048 , 
MA-083, 
MA-088 

172-173 

Right ascension and declination of velocity 

in the mean of 1950.0 coordinate 

system 

VRHA2 , VDEC2 

MA-048, 

MA-083, 

MA-088 

174-175 

Right ascension and declination of velocity 
in ecliptic mean of 1950.0 coordinate 
system 

VLOS 

MA-048 , 
MA-083, 
MA-088 

176 

Velocity - LOS angle 

LONS2 

MA-048 , 
MA-083, 
MA-088 

177 

Longitude of the vehicle in the ecliptic 
mean of 1950.0 coordinate system 

ESLOS 

MA-048 , 
MA-083, 
MA-088 

178 

Angie between Earth-Sun vector and the 
LOS 

ESLOSS 

MA-048 , 
MA-083, 
MA-088 

179 

Supplement of ESLOS 

THETA 

MA-106 

180 

Angle between x -axis and the perpendicular 
to the CSM velocity vector 

BMAG 

MA-106 

181 

Magnetic field intensity (gauss) 

LMAG 

MA-106 

182 

L-shell radius (e.r.) 

ALPH10, BETA10 , PHI 10 

MA-106 

183-185 

Orientation angles of the vehicle body 


axes wrt the UEN 



TABLE 7-1. 

EXPERIMENT SUPPORT 

DATA TAPE FORMAT 

- Continued 

Parameter name(s) 

Type of 
parameter 

Parameter 

location 

Parameter definition 

LATA, 

LONA 

MA-136 

186-187 

Geodetic latitude and longitude of the 
FOV A corner 

LATB, 

LONB 

MA-I36 

188-189 

Geodetic latitude and longitude of the 
FOV B corner 

LATP, 

LONP 

MA-136 

190-191 

Geodetic latitude and longitude of the 
principal point 

LATC, 

LONC 

MA-136 

192-193 

Geodetic latitude and longitude of the 
FOV C corner 

LAID, 

LOND 

MA-136 

194-195 

Geodetic latitude and longitude of the 
P0V D corner 

SF 


MA-136 

196 

Scale factor 

DT 


MA-136 

197 

Time interval between adjacent photographs 

ALTR 


MA-136 

198 

Attitude rate wrt the principal point 

HV 


MA-136 

199 

Horizontal velocity wrt the principal point 

TILTAZ 

MA-136 

200 

Tilt azimuth angle 

TILT 


MA-136 

201 

Tilt angle 

SELP, 

SAZP 

MA-136 

202-203 

Sun elevation and azimuth wrt the 
principal point 

LATS, 

LONS 

MA-136 

204-205 

Subsolar point - geodetic latitude and 
longitude 

ALPHA 


MA-136 

206 

Alpha angle 



TABLE 7-1. - EXPERIMENT SUPPORT DATA TAPE FORMAT - Continued 


Type of Parameter 


Parameter name(s) 

parameter 

location 

Parameter definition 

SWING 

MA-136 

207 

Swing angle 

EMISS 

MA-136 

208 

Emission angle 

PHASE 

MA-136 

209 

Phase angle 

NDA 

MA-136 

210 

North deviation angle 

XTILT 

MA-136 

211 

x-tilt angle 

YTILT 

MA-136 

212 

y-tilt angle 

HEAD 

MA-136 

213 

Heading angle 

SR 

MA-136 

214 

Slant range to the principal point 

OVR 

MA-136 

215 

Forward overlap ratio 

PHI, KAPPA, OMEGA 

MA-136 

216-218 

Angles that rotate the camera axes system 
into the local-horizontal system 

AL 

MA-136 

219 

Surface arc length between nadir and the 
principal point 

FL 

MA-136 

220 

Focal length of camera lens 

LOSX, LOSY, LOSZ 

MA-136 

221-223 

Direction cosines of the camera LOS 
vector in the geographic coordinate 
system 

GCll , GC12, GC13 

MA-136 

224-226 

Row 1 of the geographic to camera 
transformation 

LH11, LH12, LH13 

MA-136 

227-229 

Row 1 of the local horizontal to camera 


transformation 



TABLE 7-1.- EXPERIMENT SUPPORT DATA TAPE FORMAT - Concluded 


Parameter name(s) 

Type of 
parameter 

Parameter 

location 

Parameter definition 

GC21 , GC22, 

GC23 

MA-136 

230-232 

Row 2 of the geographic to camera 
transformation 

LH21 , LH22, 

LH23 

MA-136 

233-235 

Row 2 of the local horizontal to camera 
transformation 

GC31, GC32 , 

GC33 

MA-136 

236-238 

Row 3 of the geographic to camera 
transformation 

LH31, LH32, 

LH33 

MA-136 

239-241 

Row 3 of the local horizontal to camera 
transformation 



TABLE 7-II ALPHABETICAL CROSS-REFERENCE FOR ASTP EXPERIMENT SUPPORT PARAMETERS 




126 

E22 

8 

AETH 

127 

E23 

9 

AETM 

131 

E31 

10 

AETS 

132 

E32 

219 

AL 

133 

E33 

26 

ALF1 

33 

ECC1 

44 

ALF12 

51 

ECC12 

146 

ALF2 

69 

ECC3 

62 

ALF3 

87 

ECC4 

80 

ALF4 

208 

EMISS 

206 

ALPHA 

178 

ESLOS 

183 

ALPH10 

179 

ESLOSS 

109 

ALT4 

155 

F3DCA1 

198 

ALTR 

157 

F3DCB1 

29 

AZ1 

159 

F3DCC1 

47 

AZ12 

161 

F3DCD1 

65 

AZ3 

154 

F3RHA1 



156 

F3RHB1 

83 

AZ4 

158 

F3RHC1 

184 BETA10 

160 

F3RHD1 

181 

BMAG 

163 

F8DCA1 

23 

BTA1 

165 

F8DCB1 

46 

BTA12 

167 

F8DCC1 

64 

BTA3 

169 

F8DCD1 

82 

BTA4 

162 

F8RHA1 

101 

CDUX 

164 

F8RHB1 

102 

CDUY 

166 

F8RHC1 

103 

CDUZ 

168 

F8RHD1 

14 

CTEH 

220 

FL 

15 

CTEM 

224 

GCll 

16 

CTES 

225 

GC12 

27 

DLT1 

226 

GC13 

45 

DLT12 

.230 

GC21 

147 

DLT2 

231 

GC22 

63 

DLT3 

232 

GC23 

81 

DLT4 

236 

GC31 

197 

DT 

237 

GC32 

119 

EU 

238 

GC33 

120 

E12 

11 

GETH 

121 

E13 

12 

GETM 

125 

E21 

13 

GETS 


115 

GHA 

lo7 

LONA 

134 

GIMB 

189 

LONB 

3 

GMTD 

193 

LONC 

4 

GMTH 

195 

LOND 

5 

GMTM 

191 

LONP 

2 

GMTMO 

205 

LONS 

6 

GMTS 

177 

LONS2 

1 

GMTY 

221 

LOSX 

213 

HEAD 

222 

LOSY 

199 

HV 

223 

LOSZ 

34 

INCl 

210 

NDA 

52 

INC12 

35 

NODI 

70 

INC3 

53 

NOD12 

88 

INC4 

71 

N0D3 

217 

KAPPA 

89 

NOD4 

171 

L1ADC1 

218 

OMEGA 

170 

L1ARA1 

36 

OMG1 

152 

LlAZ 

54 

OMG12 

149 

L1DEC1 

72 

OMG3 

151 

L1DEC2 

90 

OMG4 

153 

LlFL 

7 

OPFLAG 

148 

L1RHA1 

215 

OVR 

150 

L1RHA2 

17 

PETH 

107 

LAT4 

18 

PETM 

186 

LATA 

19 

PETS 

188 

LATB 

209 

PHASE 

192 

LATC 

216 

PHI 

194 

LATD 

185 

PHI 10 

190 

LATP 

30 

R1 

204 

LATS 

48 

R12 

227 

LHll 

66 

R3 

228 

LH12 

84 

R4 

229 

LH13 

135 

REV 

233 

LH21 

98 

RFll 

234 

LH22 

99 

RF12 

235 

LH23 

100 

RF13 

239 

LH31 

104 

RF21 

240 

LH32 

105 

RF22 

241 

LH33 

106 

RF23 

182 

LMAG 

110 

RF31 

108 

LON4 

111 

RF32 


VRHA2 

XI 

XII 
XI 2 
X3 
X4 

XANG 
XD1 
XD12 
XD3 
XD4 
XDS1 
XS1 
XTILT 
Y1 
Yll 
Y12 
Y3 
Y4 
YANG 
YD1 • 
YD12 
YD3 
YD4 
YDS I 
YS1 
YTILT 
Z1 
Zll 
ZI2 
Z3 
Z4 
ZD1 
ZD12 
ZD3 
79 ZD4 

85 V4 97 ZDSL 

173 VDECL 94 ZSI 

175 VDEC2 

176 VLOS 
172 VRHA1 


112 RF33 174 

114 SAZIO 20 

203 SAZP 142 

145 SCSA 38 

113 SEL10 56 

202 SELP 74 

196 SF 140 

32 SMA1 23 

50 SMA12 41 

68 SMA3 59 

86 SMA4 77 

214 SR 95 

136 SUNF 92 

137 SUNH 211 

138 SUNM 21 

139 SUNS 143 

207 SWING 39 

116 Til 57 

117 T12 75 

118 T13 141 

122 T21 24 

123 T22 42 

124 T23 60 

128 T31 78 

129 T32 96 

130 T33 93 

37 TA1 212 

55 TA12 22 

73 TA3 144 

91 TA4 40 

180 THETA 58 

201 TILT 76 

200 TILTAZ 25 

31 VI 43 

49 V12 61 

67 V3 
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Figure 7-1.- Listing/microfilm format 



At the end of each run, these data will be written out on the second file or the 
ASEP tape. The data will also be listed in the following tabular format. 

Second Third Fourth 

Parameter Number Mean Moment Moment Moment 

As with the header record data, these data will appear only on the JSC master tapes. 
They will not be included on PI CCT's, 

7.4 Details of Tape Formats 

A detailed discussion of experiment support tape formats can be found in 
Appendix B. 



8.0 TEST PLAN 


The purpose of this test plan is to outline the procedures associated with 
a "closed loop" test that was conducted by JSC/TRW Task 308 in conjunction with 
the ASTP principal investigators who require ephemeris, attitude, pointing, or 
field-of-view support data, and in accordance with the Postflight Requirement Forms. 
The purpose of the test was twofold: (l) Verify that each principal investigator 

can read his computer compatible tape, and (2) verify that the experiment support 
parameters are computed correctly and displayed properly. Section 8.1 describes 
the procedures by which the test data were generated, and section 8.2 defines the 
test data formats. 


8.1 Test Data Generation Procedures 

Test data - computer compatible tapes, microfilm, and hard copies - were 
generated for each principal investigator who required ephemeris, attitude, point- 
ing, or field-of-view data. The test data were based on the ASTP premission 
reference attitude and trajectory (ref. 2). The attitude data were supplied 
by JSC/TRW Task 302, whereas the trajectory data were supplied by the Mission 
Planning and Analysis Division of JSC. 

Two segments of test data were generated. 

• Segment 1: This '5 hour and 17 minute segment (l0^ll+ m to lOP^Sl 1 ”) contains 

approximately 90 points at intervals of 0.1 to 15 minutes. 

• Segment 2: This 10 hour and 35 minute segment (l48^10 m to 158^1* 5 m ) contains 

approximately 185 points at intervals of 0.1 to 15 minutes. 

Each experiment has at least one scheduled operation period in one of these 
two intervals (see table 8-1 ). Thus, any data during a scheduled operation period 
should provide realistic values of the experiment-specific parameters for that 
experiment. For a detailed description of the attitude time line, see reference 2. 
Note that MA-089 is not listed in table 8-1 since it only requires ephemeris 
support data. 

The data generation procedure is illustrated by figure 8-1. First, the 
required attitude data were stripped from the reference attitude tape in order to 
produce a gimbal data tape compatible with the GEDIT program. (See sec. 3.1 for 
a description of the various processing functions and the associated software.) 

A parallel operation produced the trajectory tape based on the reference trajectory. 
The trajectory and attitude data were merged by the GEDIT program, to produce the 
input tape for the Apollo-Soyuz Experiment Parameter (ASEP) program, which generates 
the required experiment support data. The final step was to convert the Univac 1108 
binary tapes (ASEP output tapes) into the various computer compatible formats used 
by the principal investigators . 

8.2 Test Data Formats 

The purpose of this section is to specify the test data formats. The tape 
record format, the microfilm format, and the hard copy format are defined in 
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TABLE 8-1.- EXPERIMENT OPERATION PERIODS IN THE 
TWO SEGMENTS OF TEST DATA 


Experiment 

Operation period(s) 
in segment 1' 

Operation period(s) 
in segment 2 

MA-007 

No 

Yes 

MA-048 

No 

Yes 

MA-059 

Yes 

No 

MA-083 

Yes 

No 

MA-O88 

No 

Yes 

a MA-106 

No 

Yes 

MA-136 

Yes 

Yes 


a The MA-'10T experiment-specific parameter set is the same as the MA-106 set. 
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Figure 8-1.- Test data generation procedure. 
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detail in section 7*°> vhereas the computer compatible tape format is specified 
in table Q-II as a function of experiment. All experiments except MA-107 will 
be provided a computer compatible tape, microfilm, and a hard copy. MA-107 will 
receive only microfilm and a hard copy. 

In reading the data tapes , the investigator should design his software so that 
he can edit data by time or record number. In addition, he should be able to test 
on two important code words. These are described below. 

• QPFLAG (parameter 7) - This seven-digit positional code word specifies 
whether an experiment is operating or not. It is defined as follows. 

12 3^567 

t f t t + t t 

MA-007 MA-048 MA-059 MA-083 MA-088 MA-106 MA-136 

where a zero indicates the experiment is not operating and all leading 
zeros are dropped. For example, the number IO3OO67. specifies that MA-007, 
MA-059, MA-106, and MA-136 are operating whereas the other experiments 
are not. Likewise, the number 204500. indicates that MA-048, MA-083, and 
MA-088 are operating whereas the other experiments are not. If an 
experiment is not operating, the value of the experiment-specific parameters 
will be zero. For the test data, all experiments will be assumed to be 
operating. 

• GIMB (parameter 134) - The gimbal data quality flag denotes the status 
of the gimbal data. 

Good data 

Interpolated data - "user beware" 

No data 
Bad data 

Premission attitude reference data (MA-048 only) 

When no gimbal data are available or the data are bad, the value of all 
attitude, pointing, and FOV parameters is 0. 77777770+7* 

For the MA-136 experiment, a value of 0.88888888+8 for the geodetic location 
of an FOV corner point(s) or the principal point indicates that the associated 
line-of-sight is "above the horizon." If this is the case, then any other 
parameter that depends on the principal point or the FOV corner point(s) will be 
set equal to 0.88888888+8. 
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TABLE 8-II COMPUTER COMPATIBLE TAPE DESCRIPTION 


Experiment 

Computer 

compatible 

tape 

Tape 

density 

Tape 

parity 

Number 
of tracks 

Word 

length 

Type of 
format 

Comments 

MA-007 

(LARC) 

CDC 6600 

800 

Odd 

7 

60 

Floating 

point 

None 

MA-007 

CU. of Wyoming) 

CDC 3000 

800 

Odd 

7 

48 

Floating 

point 

JSC doe 8 not have 
Univac 1108 to 
Sigma 7 con- 
version software. 
However, Univer- 
sity of Wyoming 
does have CDC 3000 
to Sigma 7 con- 
version package. 

MA-048 

CDC 3800 

556 

Odd 

7 

48 

Floating 

point 

None 

MA-059 (JSC) 

Univac 1110 

800 

Odd 

7 

36 

Floating 

point 

None 

MA-059 

(U. of Michigan) 

IBM 360 

800 

Odd 

7 

32 

Floating 

point 

IBM 360 tape is com- 
patible with IBM 371-68 
computer . 

MA-083 

IBM 360 

556 

Odd 

7 

32 

Floating 

point 

IBM 360 tape is com- 
patible with INTER- 
DATA 70 computer. 

MA-088 

IBM 360 

556 

Odd 

7 

32 

Floating 

point 

IBM 360 tape is com- 
patible with INTER- 
DATA 70 computer. 

MA-089 

CDC 6400 

800 

Odd 

7 

60 

Floating 

point 

None 

MA-106 

CDC 6600 

800 

Odd 

7 

60 

Floating 

point 

CDC 6600 tape is 
compatible with 
CDC 7600 computer. 

MA-136 

Univac 1108 

800 

Odd 

7 

36 

Floating 

point 

Two tapes are 
required . 


9.0 TEST RESULTS 


The test data were delivered to Pi’s in two phases. First, a dummy tape was 
shipped to each PI to verily its compatibility with his computer. This tape was 
shipped March 13, 1975. These tapes were followed by the final test tapes, as 
described in section 8, on April 25, 1975. This was 1 week earlier than the PRF 
schedule. Table 9-1 summarizes results as of June 1, 1975, the test cutoff date 
as established in the PRF schedule for each experiment. Note that several Pi’s 
have not responded on one or both phases of the test. 

As of June 1, 1975, no further requests for changes in ASTP experiment support 
output can be accepted other than through the normal channel of a PRF change request, 
with review by the DEAD CCB. 
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TABLE 9-1.- VALIDATION TEST RESULTS 
Compatibility 


Experimenter 

test 

Final test 

MA-007 



LRC 

OK 

OK 

University of Wyoming 

OK 

OK 

MA-048 

OK 

No response 

MA-059 



JSC 

OK 

No response 

University of Michigan 

No response 

No response 

MA-083 

OK 

OK 

MA-088 

OK 

OK 

MA-089 

Question - 5/18/75 

No response 


No further response 


MA— 106 

OK 

No response 

MA-I07 

Did not participate 

Did not participate 

MA-I36 

OK 

OK 



10.0 QUALITY ASSURANCE FOR ASTP EXPERIMENT SUPPORT PRODUCTS 


A key element in the generation of large amounts of experiment support products 
is the quality of such data. Inaccurate and incorrect data products or untimely 
delays in the delivery of support data due to reruns can seriously jeopardize 
schedules or invalidate the scientific results. For the ASTP mission, it is 
estimated that approximately 1/2 million records of support data will be generated 
by JSC/TRW task 308 in support of principal investigators. 

The purpose of this section is to identify the controls that will be instituted 
both prior to the mission and after the mission to help guarantee high-quality data 
support products that will be generated in accordance with the production plan 
(see sec. 3.0). A summary of these controls is presented below. 


• Independent verification will be performed for all new or modified soft- 
ware prior to the mission. 

• End-to-end testing of all data processing functions will be performed 
prior to the flight with simulated data and after the mission with real 
data. 

• Independent verification will be performed for all intermediate data 
products . 


• Quality control will be maintained on final support products . 

Naturally, there are some areas that are beyond the control of the task. 
These areas include the quality and amount of input telemetry data, the inherent 
limitation in the accuracy of trajectory, the attitude and pointing input data, 
and the transmittal of the required data products. The following subsections 
detail the controls in the four areas described above. 


10.1 Independent Software Verification 

The modified preprocessing software - GREP, GEDIT, and MA-059 utility 
programs - will be independently verified prior to the mission. The Apollo-Soyuz 
Experiment Parameter (ASEP) program, a new program, will be hand checked to 
eight places. After this is successfully completed, an independent, engineering 
verification will be performed to ensure that the computation of the output 
parameters satisfies the principal investigator ( Pi) requirements. 

10.2 End-to-End-Testing 

Two types of end-to-end testing will be performed; one using simulated data 
prior to the mission and one using real data after the mission. These two tests 
are described in more detail below. 

Prior to the mission, in accordance with the DSAD PRF’s, a computer compatible 
support data tape will be delivered to each PI. This tape will be generated by the 
ASEP program with the use of ASTP reference trajectory and attitude data. Both 
nominal and experiment-specific parameters will be computed for each of the eight 
experiments supported by the task. The closed-loop test will be completed when 
each PI successfully reads the test tape and validates the test data. 
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As soon as possible after the mission, support data for a series of data takes 
over well-known U.S. sites will be generated for the MA-136 PI. This closed-loop 
test will provide the PI with a measure of the misalinement of the camera and allow 
for a final end-to-end test of the data processing system prior to the generation 
of the mission support data tapes. 


10.3 Independent Test of Intermediate Data Products 

The general operational philosophy will be to use standard data input decks, 
maintain accurate processing records, maintain a file of all intermediate data 
products, and double check data at all intermediate steps. Specifically, the 
following steps will be taken. 

Refined state vectors 

• Tracking data residual statistics from the HOPE program will be checked 
for "goodness of fit." 

• The time base will be checked for consistency with other types of data. 

• Trajectory "match point" comparisons will be made for adjacent fit 
segments as a measure of the consistency of the state vectors. 

Gimbal data 

• Gaps where gimbal data should exist will be identified. A determination 
of the status of the data will be made. 

• Segments where gimbal data exist but are not supposed to exist will be 
identified, and a determination of the validity of the data will be made. 

• Segments where the gimbal data are bad will be identified. 

• Low-bit-rate intervals will be identified. 

• Segments where the gimbal data are interpolated will be identified. 

• Plots of the first differences of the gimbal data will be made in order 

to edit "wild points." 

Master input data 

Before the generation of the support data by the ASEP program, the follow- 
ing tests will be made on the merged input data (trajectory and attitude). 

• Spot checks of the quality of the merged data will be made, 

• Data record counts will be tabulated for later comparison with the 

output record counts. This is a check on an inadvertent deletion of 
data. 
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10. 4 Quality Control on Final Support Products 

The operational philosophy will be to use standard input decks and standard 
production operations and to double check all card inputs. In addition, other 
specific measures will be taken to maintain the quality of the output products. 

• Quick-look support data for MA-136 will be generated. 

• Spot checks of the quality of the support data will be made. 

• Parameter statistics - mean and 2nd, 3rd, and ifth moments about the mean 
will be generated to help isolate problem segments. 

• Output data record count will be recorded and compared with the merged 
input record count. 

Note that verification tests will also be made on the reformatted tapes that will 
be delivered to the Pi's. 
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11.0 MISSION DESCRIPTION 


(TO BE PUBLISHED) 
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12.0 AS TP MISSION DELIVERABLES 


(TO BE PUBLISHED) 
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13.0 ARCHIVED PRODUCTS 


(TO BE PUBLISHED) 
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APPENDIX A - ASTRONOMICAL COORDINATE SYSTEMS 


1 APPENDIX A 

ASTRONOMICAL COORDINATE SYSTEMS 
A1.0 INTRODUCTION 


The purpose of this appendix is to provide an engineering description of the 
astronomical coordinate transformations used n the HOPE program. 


All transformations are performed within the framework of the celestial 
sphere. The mean-of-date coordinate system is used as a pivot system; all rota- 
tions pass through this frame even if they do not start or end there. Rotations 
are set up to go from the mean-of-date coordinate system to the input (sec. 2) 
and output (sec. 7) coordinate systems. 

Each of the various astronomical coordinate systems used in the HOPE program 
can be classified as one of the following types. 

Planetographic 
Selenographic 
Planetocentri c 
Selenocentric 
Heliocentric 

Planetographic or selenographic coordinates are coordinates that are rigidly 
attached to the planet or Moon, defined by the planet's or Moon's equator and 
prime meridian. 

Planetocentri c means referred to the center of a planet as dynamical center 
or origin of coordinates. The possible base planes for a planetocentric coordinate 
system are as follows . 

Celestial equator 
Ecliptic plane 
Planetary equator 

The possible X-axes are as follows. 

Mean equinox of date 
True equinox of date 

2 

Mean equinox of epoch 
Planetary vernal equinox 

‘‘‘The data contained in appendix A were obtained from reference b. 

2 

Epoch may be any' fixed date; for example, the beginning of the nearest 
Besselian year. 
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Selenocentric means referred to the center of the Moon. Heliocentric means 
that the center of the Sun is the origin of the coordinate system. For both these 
systems, a planetary equator is not a possible base plane. 

Section A2 of this appendix discusses the basic concepts underlying the theory 
of astronomical coordinate transformations. 

The transformation logic is discussed in section A3. 

A discussion of the computation of mean angles, precession, and physical libra™ 
tion terms and angles that define the orientation of a planetographic or planeto- 
centric coordinate system is in section AU. 


A2.0 DISCUSSION OF BASIC CONCEPTS 


A2.1 Motion of the Earth 

The equinox T is defined by the intersection of the planes of the Earth's 
equator and the ecliptic. In turn, the equator is defined as being normal to the 
Earth's pole. The primary motion of the equinox is called precession and is due 
mainly to the precession of the Earth's pole. The precessional motion of the mean 
equinox is due to the combined motions of the two planes (equator and ecliptic) 
that define it. 

The motion of the celestial pole, or of the equator, is due to the gravitational 
action of the Sun and Moon on the equatorial bulge of the Earth. It consists of 
two components, lunisolar precession and nutation. Lunisolar precession is the 
smooth long-period motion of the mean pole of the equator (the word mean indicates 
that nutation is being neglected) around the pole of the ecliptic, with an ampli- 
tude of approximately 23.5 degrees and a period of about 2 6 000 years. Nutation 
is a relatively short period motion that carries the actual (or true) pole around 
the mean pole in a somewhat irregular curve, with an amplitude of about 9 seconds 
of arc and a period of about 18.6 years. The motion of the ecliptic - that is, 
the mean plane of the Earth's orbit - is due to the gravitational action of the 
planets on the Earth and consists of a slow rotation of the ecliptic. This motion 
is known as planetary precession and gives a precession of the equinox of about 
12 seconds of arc per century and a decrease of the obliquity of the ecliptic 
(the angle between the ecliptic and the Earth’s equator) of about 47 seconds of arc 
per century. The motions of the mean pole, together with the motion of the ecliptic, 
are known as general precession. General precession is discussed further in 
section A2.1.1, and nutation is discussed in section A2.1.2. 

A2.1.1 General precession of the pole of the Earth .- As mentioned in 
section A2.1, the precession of the Earth consists of two components, lunisolar 
precession and planetary precession. It is customary to consider the two components 
together as a single quantity called general precession, which is described by 
the angles £ o , 0, Z. These precession angles are illustrated in figure A-l. The 

transformation from the mean-of-nearest-Besselian-year (mean-of-epoch) coordinate 
system to the mean-of-date coordinate system is accomplished by the [P] matrix, 
which is made up of three rotation matrices as follows. 
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e 



MOFD 

MOFNBY 


angle from mean equinox of epoch to the intersection 
of the mean equator of date and the mean equator of 
epoch 

angle between the mean equator of date and the mean 
equator of epoch 

angle from intersection of the mean equator of date 
with the mean equator of epoch to the mean equinox of 
date, measured clockwise in the mean equatorial plane 
mean-of-date coordinate system 

mean-of-nearest-Besselian-year coordinate system 


Figure A-l Geometry for general precession terms. 
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• A counterclockwise rotation about the Z-axis through the angle 


■n_ 

2 


S 


o 


c*i] - 


C0S (f - Co) 
-sin (j - C 0 ) 
0 


sin (j - C 0 ) 0 

cos (j - C 0 > 0 

0 1. 


• A counterclockwise rotation about the new X-axis through the angle 0 


[R 2 3 


10 0 
0 cos 0 sin 0 

0 -sin 9 cos 9 


• A clockwise rotation about the new Z-axis through the angle 



[R 3 1 


cos (~ + Z) 


sin + Z) 


■sin + Z) 


cos (j + Z) 


Then the [P] matrix is given by 

[P] = [r 3 ] [r 2 ] [r 1 3 

where 


Pll = cos Co cos 6 cos Z -sin Q 0 sin Z 
Pi 2 ” -sin C 0 cos 8 cos Z - cos Q q sin Z 

= -sin 9 cos Z 

P 2 i " C0S C c cos 6 sin Z + sin cos Z 


(Al) 
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p 22 

= “Sin £ o 

cos 0 

P 23 

= -sin 0 

sin Z 

P 31 

= cos r 

sin 0 

P 32 

= -sin C 

sin 9 

p 33 

= cos 0 



Thus, if r is a vector in the mean-of-nearest-Besselian-year coordinate 
system, the same vector measured in the mean-of-date coordinate system is given by 

— V 

R. 


R = [P]r 


The transpose of [P] performs the transformation from the mean-of-date coordinate 
system to the mean-of-nearest-Besselian-year coordinate system. The [P] matrix 
is orthogonal, and therefore the inverse of [P] is equal to the transpose of [?]; 

m 

that is, [P] _1 = [P] . The computation of the angles s , 0, and. Z is described 
in section kb. 2. 

A2.1.2 Nutation of the pole of the Earth .- Nutation is that part of the 
precessional motion of the pole of the Earth's equator that depends on the periodic 
motions of the Sun and Moon in their orbits around the Earth. The long-period 
motion of the mean pole is called lunisolar precession (sec. A2.l). Nutation is 
the circular motion of the true pole about the mean pole in a period of about 
19 years, with an amplitude of about 9 seconds of arc. The principal term depends 
on the longitude of the node of the Moon's orbit and has a period of 6798 days or 
18.6 years; the amplitude of this term, 9-210 seconds of arc, is known as the con- 
stant of nutation. In the complicated theory of the gravitational action of the Gun 
and Moon on the rotating nonspherical Earth, other terms arise that depend on the 
mean longitudes and mean anomalies of the Sun and Moon and on their combinations 
with the longitude of the Moon's node. The resulting shift of the true pole with 
respect to the mean pole can be resolved into corrections to the longitude (ST, nu- 
tation in longitude) and to the mean obliquity (6e, nutation in obliquity). Express- 
ions for 64* and <Se are generally developed in series form. The terms divide 
naturally into those not depending on the Moon's longitude, which can be interpolated 
at intervals of 10 days, and those depending on the Moon's longitude, with periods 
of less than about' 60 days, which cannot be so interpolated. Nutation is therefore 
divided into long-period and short-period terms. The short -period terms, which 
have periods less than 35 days, are traditionally summed separately as d¥ and 
de and are called the short-period terms of nutation in longitude and obliquity, 
respectively. The long-period terms of nutation in longitude and obliquity are 
denoted by A¥ and Ae , respectively. Thus, the total nutations in longitude 
and obliquity (fig. A-2) are given by 


APPENDIX A 


64' = A'F + d¥ 
6e = Ae + de 


The angle e in figure A-2 is the mean obliquity of the ecliptic, defined as the 
angle between the ecliptic of date and the mean celestial equator of date. The 
angle e is computed in section A4.1. The nutations and their derivatives are 
read from the JPL double-precision ephemeris tape. 

The transformation from the mean-of-date coordinate system to the true-of-date 
coordinate system is accomplished by the [N] matrix, which is made up of three 
rotation matrices as follows. 


• A counterclockwise rotation about the X-axis through an angle e 




‘i 0 0 ■ 

0 ' cos e sin e 

.0 -sin e • cos e 


(A2) 


Note: [R^ ] is the transformation from the mean-of-date coordinate system to 

the ecliptic-mean-of-date coordinate system. 

• A clockwise rotation about the new Z-axis through an angle 64' 


[R 5 3 = 


cos 
sin 6^ 
0 


-sin 6ijr 
cos 6i|f 
0 


A clockwise rotation about bhe new X-axis through an angle e + 6e 




cos (e + 6e) 
sin (e + 6e) 


0 

• sin (e + 6e) 
cos (e + 6e) 


(A3) 


Note: [Rg] is the transformation from the ecliptic-true-equinox coordinate 

system to the true-of-date coordinate system. 

Then the [n] matrix is given by 


[N] = [R 6 ][R 5 ][R 4 ] .. 


(A4) 
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where 


n ll 

— 

COS 

6f 










n 12 

=5 

-sin 

i 6^ cos 

E 









n 13 

= 

-sin 

i 6 sin 

E 









n 21 

= 

cos 

(e + 6e) 

sin 









n 22 

= 

cos 

(e + 6e) 

cos 

6li 

cos 

E 

+ sin 

(f 

-f 

6r) 

sin 

n 23 


cos 

(ft + 6 r) 

cos 

6^ 

sin 

F 

- sin 

(f 

+ 

6e) 

cos 

n 3l 


sin 

(e + 6r) 

s in 

6 + 








n 32 

as 

sin 

(ft f 6ft) 

cos 

6'|f 

cos 

F, 

- cos 

(e 

+ 

6f) 

sin 

n„„ 

33 


sin 

(f + 5e) 

cos 

6V 

sin 

E 

+ cos 

( F- 


6«) 

cos 


The transpose of [N] performs the transformation from the true-of-date coordinate 
system to the mean-of-date coordinate system. 



Figure A-2. - Geometry for nutation terms. 
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A2.2 Current Precession Calculations 

Consider an inertial mean-of-nearest-Besselian-year coordinate system that 
is coincident with the mean-of-date coordinate system; that is, the epoch date. 


where 


T t - 0.0 


T = 
ot 


^ - J£D 1900.0 
36524.219879 


(A5) 


JED = Julian ephemeris date. 

JED^qq 0 = 2415020.313, Julian ephemeris date at the beginning of the 
Besselian year 1900.0. 


Under these conditions, the precession angles C q , 0, and Z are zero at the current 
date , but their time rates of change are finite . 



= 2304.2530" 

+ 1.3973" 

T _ 
ot 

+ 0.00006" 

T 2 
ot 

(AS) 


= 2004.6850" 

- 0.8533" 

T 

ot 

- 0.00037" 

T ^ 
ot 

(aj) 

2 E 

= 2304.2530" 

+ 1.3972 " 

T 

ot 

+ 0.00006" 

T ^ 
ot 

(A8) 


The expression for the rate of change of general precession in right ascension 
nig, obtained by adding equations (a6) and (A8) , is 


niri — 


Co E + h 

♦608.506C 
tropical century) 


= 4608.5060" + 2.7945" T ot + 0.00012" T Qt 2 (seconds of arc/ 


(A9) 


The rate of change of general precession in longitude (measured in the 
ecliptic) is 


p = 0 E sin e + m^, cos e 


(A10) 
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A2.3 Greenwich Hour Angle Computation 


The angle a between the true-of-date equinox and the Greenwich meridian, 
S 

measured in the true celestial equator (fig. A-3), is given by 


a - a +6^ cos (e + 6e) (All) 

6 gm 


where 


a = angle between the mean equinox and the Greenwich meridian, measured in 

gm 

the mean celestial equator. 


The mean rotational rate of the Earth relative to the me an -of- date system 

a) is found by considering the following two components of the rotational motion. 
S 

• The daily sidereal rate or rate of Greenwich sidereal time. 

• The average daily rotation rate of seconds of arc per mean solar day. 



Figure A-3.- nutation transformation geometry. 
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The first component is found by differentiation of the Greenwich mean sidereal time 
a- 12 hours universal time, which is given by 


R ’ = l8 h 38 m 45.836 s + 8640184. 542 s T^ + 0.0929 s T* 


(A12) 


where 


T u = number of Julian centuries of 36 525 days of universal time elapsed 
since 12 hours universal time on 1900 January 0. 


The equation for T , the number of Julian centuries of 36 525 days of universal 

time elapsed from 12 hours universal time on 1900 January 0 to midnight prior to 
the current date, is 


JUD - 2415020 
T um = 36525 


(A13) 


where 


JUD = JED - — + 0.5 - 0.5, Julian universal date at midnight 

m L 86400 J 

prior to the current date 

and where 


- greatest integer function. 

AT = difference between ephemeris and universal time. 

The second computation is simply time seconds of arc/mean solar day. The sum of 
the two components is given by 



86400 + 


R 


u 

36525 


( Al4) 
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where 


R = (B6kOl8k.5k2 s ) + 0.1858 s T 
u -um 


Now, the angle a is given by 
gm 


2ir(a + si AT) 

a - 62 . 

gm 


8(d400 


( A1 ' 


where 


a 


go 


(l8 h 38 m 45.836 S + 8640184. 542 s T + 0.0929 s T 2 - I2 h ) 
' um um ■' 

R 

+ — ^ — , hours of right ascension 
36525 


A2.4 Ephemeris Time and Universal Time 

? This section is a brief discussion of the systems of time measurement that 

pertain to the determination of coordinate system transformations. 

There are three basic systems of time measurement currently in use. 

• Oidereal time - based on the diurnal motion of the stars, or the period -.if 
rotation of the Earth. 

• Universal time - based on the diurnal motion of the Sun, or the day and night 
cycle. 

• Ephemeris time - based on the orbital motion of the Earth, Moon, and Sun in 
the solar system. 

Sidereal time is directly related to the rotation of the Earth and is determined 
by observations of star transits. Equal intervals of angular motion correspond to 
equal intervals of sidereal time. Sidereal time is important because universal 
time is not determined by observation but is related to sidereal time by a numerical 
formula. Universal time, or Greenwich mean solar time, is the measure of time used 
as a basis for all civil timekeeping; it is a close approximation to the mean 
| diurnal motion of the Sun. 

Universal time is defined as 12 hours plus the Greenwich hour angle of a point 
on the equator whose right ascension, measured from the mean equinox of date, is 
given by equation (A12). Note that HOPE assumes universal time. 

Equation (A12) is used to relate rigorously sidereal and universal time. 

Thus, the ‘ sidereal and universal time systems are not independent, and the use 
i of one instead of the other is purely a matter of convenience. 
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Ephemeris time is a theoretically uniform time scale that came into being when 
it was discovered that the Earth's rotation rate is subject to slight variations. 
The length of the ephemeris second is fixed by definition. Ephemeris time is con- 
sidered to be the independent time argument of the ephemerides of the Sun, Moon, 
and planets. There is a fundamental difference between ephemeris time and the two 
systems of time measurement discussed previously. Ephemeris time is defined to 
be uniform and is determined, in practice, through the motion of the Moon in its 
orbit around the Earth. Sidereal time (or universal time) is determined from the 
actual rotation of the Earth. Thus, it is not possible to express one system 
in terms of the other; the relation between them i&ust be determined empirically. 

In the determination of coordinate transformations, ephemeris time is the 
independent time argument in the precession, nutation, and libration terms. 
Universal time is used in the computation of R u (eq. A12) . Sidereal time is not 

used as a time argument in any of the expressions used to determine coordinate 
transformations . 


A3.0 TRANSFORMATION DESCRIPTION 


This section describes the logic and philosophy for the astronomical coordinate 
transformations (rotations, not translations). The underlying philosophy is that 
the mean-of-date coordinate system is used as the pivot coordinate system; all 
rotations pass through this frame even if they do not start or end there. Rotations 
are set up to go from the mean-of-date coordinate system to the input and output 
coordinate systems. The astronomical coordinate systems are listed below and 
described in table A-I . 

The HOPE program is capable of generating the matrix for rotation between 
any two of the following sets of astronomical coordinate systems. 

(1) Mean of date 

(2) True of date 

(3) Ecliptic, true equinox 

(4) True sele nographic 

( 5 ) Mean of epoch (where "epoch" may be any fixed date; for example, the 
beginning of the nearest Besselian year) 

(6) Mean selenographic 

( 7 ) Ecliptic, mean equinox 

(8) Planetographic 

(9) Planetocentric, planetary equator 

In the coordinate systems listed above, the vector is represented in 
Cartesian form (table A-l). 
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TABLE A-I CARTESIAN COORDINATE SYSTEMS 


Number 

Name 

Base plane 

X-axis 

(1) 

Mean of date 

Mean celestial equator 
of date 

Mean equinox of date 

(2) 

True of date 

True celestial equator of 
date 

True equinox of date 

(3) 

Ecliptic , true 
equinox 

Ecliptic plane of date 

True equinox of date 

(4) 

True selenographic 

True lunar equator 

Prime lunar meridian 

(5) 

Mean of epoch 

Mean celestial equator of 
epoch 

Mean equinox of epoch 

(6) 

Mean selenographic 

Mean lunar equator of date 

Prime lunar meridian 

(7) 

Ecliptic, mean 
equinox 

Ecliptic plane of date 

Mean equinox of date 

(8) 

Planetographi c 

Planet's equator of date 

Planet's prime meridian 

(9) 

Planetocentric , 
planetary equator 

Planet's equator of date 

Planetary vernal 
equinox of date a 


a De fined at this time to he the mean planetary vernal equinox of date. 

By noting the following definitions - 

e = mean obliquity. 

6V = total nutation in longitude. 

6e = total nutation in obliquity, 
ft = longitude of the mean ascending node of the lunar orbit, 
a = physical libration in node. 

I = mean inclination of the Moon's equatorial plane with respect to the 
ecliptic plane. 

p = physical libration in inclination. 

(J = geocentric mean longitude of the Moon. 
t = physical libration in longitude. 

Z + — = angle from intersection of the mean e qua .or of date with the mean 

equator of epoch to the mean equinox of date, measured clockwise in 
the mean equatorial plane. 

6 = angle between the mean equator of date and the mean equator of epoch. 

~ angle from mean equinox of epoch to the intersection of the mean equator 
_ , of datg^ and the, mean equator of epoch . 
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a = right ascension of the planet's north pole on the celestial sphere, 

^ referred to the mean pole and equinox of date. 

6 = declination of the planet's north pole on the celestial sphere, referred 

^ to the mean pole and equinox of date. 

V = hour angle of the planetary vernal equinox, measured from the planetary 
^ prime meridian. 

4> = angle between the planetary vernal equinox and the ascending nocr-„* of 
the planetary orbit, measured counterclockwise along the planetary 
equator. 

and using coordinate system (l) as the pivot frame, the various transformations are 


summarized as 

follows . 




Option 

Rotation 1 
about X 

Rotation 2 
about Z 

Rotation 3 
about X 

Rotation 4 
ab out Z 

(1) ■* (2) 

e 

-6¥ 

-( £+<$£ ) 

0 

(1) + (3) 

e 

-6V 

0 

0 

(1) + (4) 

e 

£2 + a + it 

I + p 

(C-ft) + (t-i 

(1) - (5) 

0 

z + f 

-0 

S 2 

(1) +■(6) 

£ 

ft + ir 

I 

C - ft 

(1) + (7) 

£ 

0 

0 

0 

(1) - (8) 

0 

, TT 

“p + 2 

J - 6 

2 p 

v - * 

p 

(1) (9) 

0 

a + ~ 

p 2 

7r JC 

2 ~ 6 P 



The logic for the transformations between coordinate systems is shown in 
figure A- 4. 
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The arrows between coordinate systems indicate the procedure by which 
vectors are transformed. For example, a transformation between coor- 
dinate systems (5) and (6) is not made directly; the procedure is 
(5) (l) (6) or vice versa. 

Figure A-k.- Logic diagram for astronomical coordinate 
system transformations. 


Al*.0 BASIC COMPUTATIONS 


Ah.l Computation of Mean Angles 

There are eight mean angles to be determined: six mean longitudes, the mean 

obliquity, and the angle of inclination of planetary orbit to ecliptic. Hie mean 
longitudes , the mean obliquity, and the angle of inclination of planetary orbit 
plane to ecliptic are defined as follows. 

f = mean longitude of the Sun's perigee measured in the ecliptic plane from 
the mean equinox of date. 

f = mean longitude of the Moon's perigee measured in the ecliptic plane from 
the mean equinox of date to the mean ascending node of the lunar orbit 
and then along the orbit. 

L = mean longitude of the Sun measured in the ecliptic plane from .the menu 
equinox of date. 

ft = longitude of the mean ascending node of the lunar orbit measured in the 
ecliptic plane from the mean equinox of date * 
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C - geocentric mean longitude of the Moon measured in the ecliptic plane 
from the mean equinox of date to the mean ascending node of the lunar 
orbit and then along the orbit. 

ft = angle, measured along the ecliptic, between the mean (Earth) vernal 

^ equinox and the ascending node of the planetary orbit on the plane of 

the ecliptic, 

= inclination of the planetary orbit plane to the ecliptic. 

e = angle between the mean celestial equator and the ecliptic (mean obliquity) 

The mean longitude of the Sun's perigee, T, and the mean obliquity, e, are 
slowly varying functions of time as compared to T' t L, ft, and The nutations 

and librations are small periodic variations about the mean position. 

The expressions for the mean angles and the mean obliquity are obtained 
through the use of the following algorithm. 

T„ = - - D -- ~ — 2* 1 (A1 6) 


where 

Tg = time in Julian centuries of 3 6 525.0 ephemeris days from the epoch date. 
J.D. = Julian date, ephemeris time. 

and where 

epoch date = (J.D. ) 1900 . 0 Jamiaiy 0>5 = 2415020.0 

The mean angles and the mean obliquity are then given by the following 


expressions . 

T = 1012395.00" + 6189.03" T_ + 1.63" T|, + 0.012" T^ (A17) 

hi hi 

r' = 1203586.40 " + 14648522.52" t £ - 37.17 " t| - 0.045 " t| (ai8) 

L = 1006908.04" + 129602768.13" T_ + 1.089 " t2 (A19) 

& ill 

ft = 933059.79" - 6962911.23" T w + 7.48" t| + 0.008" (A20) 

C = 973562.99 " + 1732564379.31" t £ - 4.08" t| + 0.0068 " t| (A 21) 
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£2 = 175631.19" + 2775.57" T_ - 0.005" - 0.0192" T l (A22) 

p E E E 

e = 84428.26” - 46.845" T_ - 0,0059" t£ + 0.00181" (A23) 

ill 

e = 6661.20” - 2.430" T_ + 0.0454" T* (A24) 

p EE 

A4.2 Computation of Precession Terms 

This section presents expressions for the determination of the three terms 
? , Z, and 0. These three terms determine the general precession of the equinox 

between two dates . 

Each of the equations for the components of the total precession are expressed 
in terras of the following quantities . 

'T ^ = time from 1900.0 to the selected reference epoch in tropical centuries. 
Tj. = elapsed time since the reference epoch in tropical centuries. 

The quantity T is a constant determined by the particular reference epoch 

chosen. The quantity T is computed according to equation (A5). 

ot 

The quantity T^ is computed according to 

rn JED - epoch 

t = 36524.219879 

The three general precession terms are defined as follows. 

C Q + Z = general precession in right ascension. 

0 = precession in declination. 

7T 

- X, - angle from mean equinox of epoch to ascending node of mean equator 
0 of date. 

The values of the above variables are determined from the following equations . 
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(2304.2530" + 1.3973" T . + 0.00006" T 2 ) T 
V ot ot / t 

+ ( 0.3023” - 0.00027" T J T 2 + 0.0180" T 3 
\ Ot' t t 


(2004.6850 
- (0.4267 

(2304.2530 
+ (l.0950 


» - 0.8533" T - 0.00037" T } ) T 
ot ot t 

" + 0.00037" T ) T 2 - 0.0418" T 3 
ot/ t t 


" + 1.3972" T + 0.00006" T 2 ) T. 

ot ot ' t 

» + 0.00039" T \ T 2 + 0.01832" T ' 
ot/ t t 


(A25) 


(A26) 


(A27) 




APPENDIX B - EXPERIMENT SUPPORT TAPE FORMATS 
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EXPERIMENT SUPPORT TAPE FORMATS 
B1.0 INTRODUCTION 


The objective of this appendix is to provide a summary of the formats of the 
tapes associated with ASTP experiment support data processing. The tape formats 
that are summarized in the following sections include: 

• Computer word tape 

• Camera shutter tape 

• ASEP input tape 

• ASEP output tape 

• MA-059 observation tape 

• MA-059/HOPE output tape 

• MA-059 experiment support tape 

All tapes are 7-track, 800 bpi Univac 1108 compatible. 


B2.0 COMPUTER WORD TAPE FORMAT 


The computer word tape contains the downlink telemetry parameters that are 
required for the generation of ASTP experiment support data. The computer word 
tape is an input to the GEDIT program (see sec. 2.2.1). The tape has no header 
records; the data records consist of 72 single-precision words in engineering units. 
Table B-I specifies the parameter name, the location on tape for each parameter, 
and the computer word number associated with the same parameter on the Skylab computer 
word tape. 

B3.0 CAMERA SHUTTER TAPE FORMAT 

The camera shutter tape contains the times of the "shutter open" event for the 
two MA-136 cameras. The tape contains two header records; each data record contains 
60 integer words. 


PRECEDING EEOE 


s® HL ““ 
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TABLE 

B-I . 

- PARAMETER 

LIST FOR 

THE IDSD COMPUTER 

WORD TAPE 


Output 

parameter 

location 

Computer 
word no. 

Parameter 


Output 

parameter 

location 

Computer 
word no . 

Parameter 


Output 

parameter 

location 

Computer 
word no. 

I 

1 

RANGE TIME 


25 

98 

MARICDUZ 


49 

177 

2 

26 

DELV 

X 

26 

99 

MAROCDUT 


50 

178 

3 

27 

DELV 

Y 

27 

100 

MARICDUX 


51 

200 

4 

28 

DELV 

Z 

28 

101 

CDUT 


52 

207 

5 

32 

RCSM 

X 

29 

102 

CDUS 


53 

208 

6 

33 

RC»t 

Y 

30 

104 

REFMAT 

0 

54 

209 

7 

34 

RCJM 

Z 

31 

105 

REFMAT 

1 

55 

210 

8 

35 

VC91 

X 

32 

106 

REEMAT 

2 

56 

211 

9 

36 

VC 91 

Y 

33 

107 

REFMAT 

3 

57 

212 

10 

37 

VC 91 

Z 

34 

108 

REFMAT 

4 

58 

213 

11 

38 

T STATE 


35 

109 

REFMAT 

5 

59 

215 

12 

41 

CDU 

X 

36 

110 

REFMAT 

6 

60 

216 

13 

42 

CDU 

Y 

37 

111 

REFMAT 

7 

61 

217 

14 

43 

CDU 

Z 

38 

112 

REFMAT 

8 

62 

224 

15 

80 

R OTHER 

X 

39 

125 

ADOT 

X 

63 

238 

16 

81 

R OTHER 

Y 

40 

126 

ADOT 

Y 

64 

239 

17 

82 

R OTHER 

Z 

41 

127 

ADOT 

Z 

65 

240 

18 

83 

V OTHER 

X 

42 

156 

TNOW 


66 

247 

19 

84 

V OTHER 

Y 

43 

170 

PI PA 

X 

67 

248 

20 

85 

V OTHER 

Z 

44 

171 

P1PA 

Y 

68 

259 

21 

86 

T OTHER 


45 

172 

PI PA 

Z 

69 

269 

22 

95 

MARKDOWN 


46 

174 

CDU 

X 

70 

270 

23 

96' 

MARICDUY 


47 

175 

CDU 

Y 

71 

508 

24 

97 

MAROCDUS 


48 

176 

CDU 

Z 

72 

509 


Parameter 
VHF CUT 
VHF TIME 
T PIP 
MARK2TIME 
MARK2 ICDUY 
MARK2PCDUS 
MARK2ICDUZ 
MARK20CDUT 
MARK2ICDUX 
CDUT 
ADOT X 
ADOT Y 
ADOT Z 
Ri 
OGC 
IGC 
MGC 
CODE 
CODE 

TRKMKCNT 
ST ARID 1 
STARID2 
DVTOTAL 
MARKTIME 
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B3.1 Header Records 

The two header records are identical and contain 48 BCD characters in the 
foil owing format . 

APOLLObSOYUZbTESTbPROJ. bEXPbMA136bbbMMDDYYbRUN bX 
MM/DD/YY = date of run, where 
MM = month (01-12) 

DD = day (01-31) 

YY = year (75) 

= run number 
b - blank 


B3.2 Data Records 

Each physical record contains 30 data frames. Each data frame contains two 
words of data, as follows. 


Word no. 


Contents 


1 Time (integer milliseconds) 

2 Shutter open data 


a. Bits 1-24: Unused (bit 1 - leftmost or MSB) 


b. Bits 25-36: CK1043 (shutter open value) 


Time that is used is the G.m.t. of the frame in which "shutter open" status occurred 
or that G.m.t. plus 10 milliseconds; which one depends on whether the "shutter open" 
status occurred with the first or second sample in the frame. The remaining words 
of an incomplete record are filled with 7' s . If the final data frame completes a 
record, then the words in the last record are filled with 7* s . 


B4.0 ASEP INPUT TAPE FORMAT 


The ASEP input tape that is produced by the GEDIT program (see sec. 2.2.1) 
provides time-ordered ephemeris, attitude, and transformation data from which all 
experiment support parameters can be computed. Table B-II specifies the format of 
this tape. 
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TABLE B-II ASEP INPUT DATA TAPE FORMAT 
(a) Initial record 


Word 

ID 

Type 


Description 

1 

BTIME 

DP 

Yr (mod. 1900) 


2 

BTIME 

DP 

Month 


3 

BTIME 

DP 

Day 

• Base date 

4 

BTIME 

DP 

Hr 


5 

BTIME 

DP 

Min 


6 

BTIME 

DP 

Sec 


7 

GET 

DP 

Min from launch 

of base time 


8 

GHAE 

DP 

Right ascension 

of Greenwich midnight day of epoch 

9 

GHAL 

DP 

Right ascension 

of Greenwich midnight day of launch 

10-20 

TITLE 

HOL 

66 alphanumeric 

characters from HOPE title card 

21-23 

N/A 

SP 

Not used 





(b) Data 

record (units: Ion, sec, and deg) 

Word 

ID 

Type 

Description 

1 

TAET 

DP 

Time (min) from BTIME 

2-4 

PECL 

DP 

Vehicle position - ecliptic mean of 1950.0 

5-7 

VECL 

DP 

Vehicle velocity - ecliptic mean of 1950.0 

8- 10 

AECL 

DP 

Vehicle acceleration - ecliptic mean of 1950.0 

11-13 

PTOD 

DP 

Vehicle position - true of date 

14-16 

VTOD 

DP 

Vehicle velocity - true of date 

17-19 

ATOD 

DP 

Vehicle acceleration - true of date 

20-22 

PGEO 

DP 

Vehicle position - geographic inertial 

23-25 

VGEO 

DP 

Vehicle velocity - geographic inertial 

26-28 

AGEO 

DP 

Vehicle acceleration - geographic inertial 

29-31 

1M50 

DP 

Vehicle position - mean of 1950.0 

32-34 

VM50 

DP 

Vehicle velocity - mean of 1950.0 

35-37 

AM50 

DP 

Vehicle acceleration - mean of 1950.0 

38-40 

PROT 

DP 

Vehicle position - geographic rotating 

41-43 

VROT 

DP 

Vehicle velocity - geographic rotating 

44-46 

AROT 

DP 

Vehicle acceleration - geographic rotating 

47-55 

EM AT 

DP 

Mean of 1950.0 to ecliptic mean of 1950.0 matrix 

56-64 

TO AT 

DP 

Mean of 1950.0 to true-of-date matrix 

65-73 

WTO AT 

DP 

Mean of 1950.0 to geographic inertial matrix 

74-76 

SUNP 

DP 

Sun position - mean of 1950,0 

77-79 

SUN V 

DP 

Sun velocity - mean of 1950.0 

80-82 

CDUX, 

CDUY , 
CDUZ 

DP 

Gimbal angles 

83-85 

MUNP 

DP 

Moon position — mean of 1950.0 

86-88 

MUNV 

DP 

Moon velocity - mean of 1950.0 

89 

GIMF 

DP 

Gimbal status flag 

90-91 

N/A 

DP 

Not used 




(c) Terminating record 

Word 



Type Value 

1-91 



DP 0.00 

92 



SP EOFAAA (where A = blank) 
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B5.0 ASEP OUTPUT TAPE FORMAT 


The ASEP output data tape consists of data records with 241 single-precision 
words per record. Each data tape is terminated with an end-of-file (EOF). The 
format for the data records is specified in table B-III. The table lists the 
parameter name, the type of parameter, the parameter location, and the parameter 
definition. Note that the parameter units are kilometers, seconds, and degrees 
unless otherwise specified. 


B6.0 MA-059 OBSERVATION TAPE FORMAT 


A special utility program is used to generate an angle-tracking data tape in 
the LOVE EDITED MASTER format based on the ASTP gimbal angle data. The format is 
described in table B-IV. 


B7.0 MA-059/H0PE OUTPUT TAPE 


This HOPE output tape contains special MA-059 parameters derived from the 
reconstruction of the Soyuz trajectory with respect to the CSM. The Fortran, 
binary, double-precision tape has three word records. There are no header records. 

Word Word description 


1 Time (minutes from AET) 


2 


Range (kilometers) 


3 Y (degrees) 

A double EOF will be placed at the end of the data records. 

B8.0 MA-059 EXPERIMENT SUPPORT TAPE 


This tape is created by merging the nominal experiment support parameters from 
the ASEP output tape with the parameters from the MA-059/HOPE output tape. The 
Fortran, binary tape has records with 139 single-precision words. There are no 
header records. 

Word Word description 


1-135 

Nominal ASEP parameters 

136 

Range (word 2 of A) 

137 

6 (word 180 of ASEP) 

138 

Y (word 3 of A) 

139 

A - 0 - Y (computed) 


A double EOF will be placed at the end of the data records. 
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Type of 

Parameter 


Parameter name(s) 

parameter 

location 

Parameter definition 

GMT 

Nominal 

1-6 

Greenwich mean time (yr (mod. 1900), month, day, hr, 
min, sec) 

OPFLAG 

Nominal 

7 

Code word to designate which experiments are operating 

AETH, AETM, AETS 

Nominal 

8-10 

Ground elapsed time (hr, min, sec) 

GETH, GETM, GETS 

Nominal 

11-13 

ASTP mission elapsed time (hr, min, sec) 

CTEH, CTEM, CTES 

Nominal 

14-16 

Onboard clock elapsed time (hr, min, sec) 

PETH, PETM, PETS 

Nominal 

17-19 

Phase elapsed time (hr, min, sec) 

XI, Yl, Z1..XD1, YDl , 
ZD1 

Nominal 

20-25 

State vector - Cartesian elements in the mean of 1950.0 
coordinate system 

ALF1 , DLTI , BTAl , AZ1, 
Rl, VI 

Nominal 

26-31 

State vector - ADBARV elements in the mean of 1950.0 
coordinate system 

SMA1, ECC1, INC1 , NODI, 
OMG1 , TAl 

Nominal 

32-37 

State vector - the Keplerian elements in the mean of 1950.0 
coordinate system 

X12 , YI2, Z12 , XD12, 
YD12, ZD12 

Nominal 

38-43 

State vector - Cartesian elements in the geographic rotating 
coordinate system 

ALF12, DLT12, BTAl 2 , 
AZ12, R12, V12 

Nominal 

44-49 

State vector - ADBARV elements in the geographic rotating 
coordinate system 

SMA12, ECC12, INC12, 
NODI 2, OMG12, TA12 

Nom inal 

50-55 

State vector - Keplerian elements in the geographic rotating 
coordinate system 

X3, Y3, Z3 , XD3 , ZD3 

Nominal 

56-61 

State vector - Cartesian elements in the true of date 
coordinate system 

ALF3 , DLT3, BTA3, AZ3, 
R3 , V3 

Nominal 

62-67 

State vector - ADBARV elements in the true of date 
coordinate system 





TABLE B-III.- EXPERIMENT SUPPORT DATA TAPE FORMAT - Continued 


Type of Parameter 


Parameter name(s) 


parameter 

location 

Parameter definition 

SMA3, ECC3 , INC3 , 
0MG3 , TA3 

NOD3, 

Nominal 

68-73 

State vector - Keplerian elements in the true of date 
coordinate system 

X4 , Y4, Z4, XD4, 
ZD4 

YD4, 

Nominal 

74-79 

State vector - Cartesian elements in the geographic inertial 
coordinate system 

ALF4, DLT4, BTA4, 
R4, V4 

AZ4, 

Nominal 

80-85 

State vector - ADBARV elements in the geographic inertial 
coordinate system 

SMA4, ECC4, INC4, 
0MG4, TA4 

NOD4, 

Nominal 

86-91 

State vector - Keplerian elements in the geographic 
coordinate system 

XS1, YS1, ZS1, XDS1, 
YDS1, ZDS1 

Nominal 

92-97 

Solar Cartesian elements in the mean of 1950.0 coordinate 
system 

RF11, RF12, RF13 


Nominal 

98-100 

Row 1 of REFSMMAT 

CDUX, CDUY, CDUZ 


Nominal 

101-103 

Gimbal angles 

RF21 , RF22, RF23 


Nominal 

104-106 

Row 2 of REFSMMAT 

LAT4, LON4, ALT4 


Nominal 

107-109 

Geodetic altitude, latitude, and longitude 

RF31 , RF32, RF33 


Nominal 

110-112 

Row 3 REFSMMAT 

SEL10, SAZ10 


Nom inal 

113-114 

Solar azimuth and elevation wrt the subvehicle point 

GHA 


Nominal 

115 

Greenwich hour angle 

Til, T12 , T13 


Nominal 

116-118 

Row 1 of the T matrix (mean of 1950.0 to true of date) 

Ell, E12, E13 


Nominal 

119-121 

Row 1 of the E matrix (mean of 1950.0 to ecliptic 
mean of 1950.0) 

T21, T22, T23 


Nominal 

122-124 

Row 2 of the T matrix 


TABLE B-III EXPERIMENT SUPPORT DATA TAPE FORMAT - Continued 


Pamameter name(s) 

Type of 
parameter 

Parameter 

location 

E21, E22, E23 

Nominal 

125-127 

T31 , T32 , T33 

Nominal 

128-130 

E31 , E32, E33 

Nominal 

131-133 

GIMB 

Nominal 

134 

REV 

Nominal 

135 

SUNF 

MA-007 

136 

SUNH, SUNM, SUNS 

MA-007 

137-139 

XANG , YANG 

MA-007 

140-141 

XU, YU, Zll 

MA-048, 
MA-083 , 
MA-088 

142-144 

SCSA 

MA-048 , 
MA-083, 
MA-088 

145 

ALF2, DLT2 

MA-048 , 
MA-083, 
MA-088 

146-147 

L1RHA1, L1DEC1 

MA-048 , 
MA-083, 
MA-088 

148-149 

L1RHA2, L1DEC2 

MA-048 , 
MA-083, 
MA-088 

150-151 


Parameter definition 
Row 2 of the E matrix 
Row 3 of the T matrix 
Row 3 of the E matrix 
Gimbal status flag 
Rev . number 
Sunrise/sunset flag 
Sunrise/sunset time (hr, min, sec) 
x-angle and y-angle deviations 
Geomagnetic position of the vehicle 

Vehicle - Sun angle 


Right ascension and declination of the vehicle in the 
ecliptic mean of 1950.0 coordinate system 


Right ascension and declination of the LOS in the mean 
of 1950.0 coordinate system 


Right ascension and declination of the LOS in the ecliptic 
mean of 1950.0 coordiante system 


TABLE B-III EXPERIMENT SUPPORT DATA TAPE FORMAT - Continued 


Parameter name(s) 

Type of 
parameter 

Parameter 

location 

L1AZ, L1EL 

MA-048 , 
MA-083, 
MA-088 

152-153 

F3RHA1 , F3DCA1 

MA-048 , 
MA-083, 
MA-088 

154-155 

F3RHB1, F3DCB1 

MA-048 , 
MA-083, 
MA-088 

156-157 

F3RHC1, F3DCC1 

MA-048 , 
MA-083, 
MA-088 

158-159 

F3RHD1, F3DCD1 

MA-048 , 
MA-083, 
MA-088 

160-161 

F8RHA1, F8DCA1 

MA-048 , 
MA-083, 
MA-088 

162-163 

F8RHB1, F8DCB1 

MA-048 , 
MA-083, 
MA-088 

164-155 

F8RHC1, F8DCC1 

MA-048 , 
MA-083, 
MA-088 

166-167 

F8RHD1, F8DCD1 

MA-048 , 
MA-083, 
MA-088 

168-169 


Parameter definition 
Azimuth and elevation of the LOS 


Right ascension and declination of FOV A in the mean 
of 1950.0 coordinate system for MA-083 


Right ascension and declination of FOV B in the mean of 

1950.0 coordinate system for MA-083 


Right ascension and declination of i-'/V C in the mean of 

1950.0 coordinate system for MA-083 


Right ascension and declination of FOV D in the mean 
of 1950.0 coordinate system for MA-083 


Right ascension and declination of FOV A in the mean of 

1950.0 coordinate system for MA-088 


Right ascension and declination of FOV B in the mean of 

1950.0 coordinate system for MA-088 


Right ascension and declination of FOV C in the mean of 

1950.0 coordinate system for MA-088 


Right ascension and declination of FOV D in the mean of 

1950.0 coordinate system for MA-088 


TABLE B-1II EXPERIMENT SUPPORT BATA TAPE FORMAT - Continued 


Parameter name(s) 

Type of 
parameter 


Parameter 

location 

iParameter definition 


L1ARA1 

, L1ADC1 

MA-048 

MA-083 

MA-088 


170-171 

Corrected right ascension and declination of 
LOS in the mean of 1950.0 coordinate system 

the 

VRHA1, 

VDECl 

MA-048 

MA-083 

MA-088 


172-173 

Right ascension and declination of velocity 
mean of 1950.0 coordinate system 

in the 

VRHA2, 

VDEC2 

MA-048 

MA-083 

MA-088 


174-175 

Right ascension and declination of velocity 
mean of 1950.0 coordinate system 

in ecliptic 

VLOS 


MA-048 

MA-083 

MA-088 


176 

Velocity - LOS angle 


LONS2 


MA-048 

MA-083 

MA-088 


177 

Longitude of the vehicle in the ecliptic mean of-, 
1950.0 coordinate system 

ESLOS 


MA-048 

MA-083 

MA-088 


178 

Angle between Earth-Sun vector and the L;v'3 


ESLOSS 


MA-048 

MA-083 

MA-088 


179 

Supplement of ESLOS 


THETA 


MA-106 


180 

Angle between X v ~axis and the perpendicular 
CSM velocity vector 

to the 

BMAG 


MA-106 


181 

Magnetic field intensity (gauss) 


LMAG 


MA-106 


182 

L-shell radius (e.r.) 
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TABLE B— III EXPERIMENT SUPPORT DATA TAPE FORMAT - Continued 


Parameter name(s) 

Type of 
parameter 

Parameter 

location 

ALPH10 , BETA10, PHI10 

MA— 106 

183-185 

LATA, LONA 

MA-136 

186-187 

LATB, LONB 

MA-136 

188-189 

LATP, LONP 

MA-136 

190-191 

LATC, LONC 

MA-136 

192-193 

LATD, LOND 

MA-136 

194-195 

SF 

MA-136 

196 

DT 

MA-136 

197 

ALTR 

MA-136 

198 

HV 

MA-136 

199 

TILTAZ 

MA-136 

200 

TILT 

MA-136 

201 

SELF, SAZP 

MA-136 

202-203 

LATS, LONS 

MA-136 

204-205 

ALPHA 

MA-136 

206 

SWING 

MA-136 

207 

EMISS 

MA-136 

208 


Parameter definition 

Orientation angled of the vehicle body axes wrt 
the UEN 

Geodetic latitude and longitude of the FOV A corner 

Geodetic latitude and longitude of the FOV B corner 

Geodetic latitude and longitude of the principal point 
Geodetic latitude and longitude of the FOV C corner 

Geodetic latitude and longitude of the FOV D corner 

Scale factor 

Time interval between adjacent photographs 
Attitude rate wrt the principal point 
Horizontal velocity wrt the principal point 
Tilt azimuth angle 
Tilt angle 

Sun elevation and azimuth wrt the principal point 
Subsolar point - geodetic latitude and longitude 
Alpha angle 
Swing angle 
Emission angle 
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TABLE B-III EXPERIMENT SUPPORT DATA TAPE FORMAT - Concluded 


Type of Parameter 


Parameter name(s) 

parameter 

location 


Parameter definition 

PHASE 

MA-136 

209 

Phase 

angle 

NDA 

MA-136 

210 

North 

deviation angle 

XTILT 

MA-136 

211 

x-tilt 

angle 

YTILT 

MA-136 

212 

y-tilt 

angle 

HEAD 

MA-136 

213 

Heading angle 

SR 

MA-136 

214 

Slant' 

range to the principal point 

OVR 

MA-136 

215 

Forward overlap ratio 


PHI, KAPPA, OMEGA 

MA-136 

216-218 

Angles that rotate the camera axes system into the 
local-horizontal system 

AL 

MA-136 

219 

Surface arc 

length between nadir 

and the principal point 

FL 

MA-136 

220 

Focal 

length 

of camera lens 


LOSX, LOSY, LOSZ 

MA-136 

221-223 

Direction cosines of the camera LOS vector in the 
geographic coordinate system 

GC11, GC12, GC13 

MA-136 

224-226 

Row 1 

of 

the 

geographic to camera 

transformation 

LHllj LH12, LH13 

MA-136 

227-229 

Row 1 

of 

the 

local horizontal to 

camera transformation 

GC21, GC22 , GC23 

MA-136 

230-232 

Row 2 

of 

the 

geographic to camera 

transformation 

LH21 , LH22 , LH23 

MA-136 

233-235 

Row 2 

of 

the 

local horizontal to 

camera transformation 

GC31, GC32, GC33 

MA-136 

236-238 

Row 3 

of 

the 

geographic to camera 

transformation 

LH31, LH32, LH33 

MA-136 

239-241 

Row 3 

of 

the 

local horizontal to 

camera transformation 
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TABLE B-IV . - MA-059 OBSERVATION TAPE FORMAT 
(a) Header record 

The header record is defined as follows. 


Word 


Value 

1-2 

Double- 

precision Julian date ( inpu?J 

3 

Date of 

run: DAY-MONTH-YEAR (12 bits . 

4 

Time of 

run: HRS-MIN-SEC (12 bits eac! 

5 

1108 

M 

6 

ASTER 


7-50 

Comments 


(b) 

Data record 

Word 

Bits 


1 

0-17 

Zero 

1 

18-35 

Sequence number (1,3,5,**.) 

2 

0-17 

CSM (alphanumeric) 

2 

18-35 

SOY (alphanumeric) 

3 

0-5 

6 (right adjusted) 

3 

6-11 

Zero 

3 

12-17 

Zero 

3 

18-23 

Zero 

3 

24-29 

2 (right adjusted) 

3 

30-35 

6 (right adjusted) 

4 

0-17 

Same as word 3, bits (0-17) 

4 

18-35 

Blank 

5 

0-35 

Bit 32 * 1; 0 otherwise 

6 

0-35 

Zero 
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TABLE B-IV . - MA-059 OBSERVATION. TAPE FORMAT - Concluded 

(b) Data record - Concluded 


Word 

Bits 


7-8 

0-71 

Time from base day 

9-10 

0-71 

Zero 

11-12 

0-71 

Zero 

13-14 

0-71 

t/ 2 rad 

15-16 

0-71 

Zero 

17-18 

0-71 

Zero 

19-20 

0-71 

Inner GA (rad) 

21-22 

0-71 

Middle GA (rad) 

23-24 

0-71 

Outer GA (rad) 

25-32 

_ 

Zero 


The total size of the data record is 640 words (twenty 32-word frames). 


(c) Final record 

The tape is terminated by an end-of-tape frame. This 32-word frame consists 
of 30 alphanumeric words of THE END in words 1-6 and 9-32 and 1 double-precision 
word (7-8) of 1.D30. This end-of-tape frame is placed in the remaining frames of 
an incomplete data record or in all 20 frames of the data record following a com- 
plete data record. 

(d) EOF 

An EOF mark is placed after the final record. 
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APPENDIX C 


COMPUTER COMPATIBLE TAPE FORMATS 


1.0 INTRODUCTION 


The purpose of this appendix is to summarize the tape formats associated with 
the conversion of ASEP output tapes from the Univac ll80 computer to other computers 
For completeness, the formats for all principal investigator experiment support 
tapes are included below. 


Machine 
Density- 
Parity 

Number tracks 
Word length 
Word type 

Number of words - 241 data words per record. See appendix B5.0 for the 
definition of the data words . 

Read - Use Fortran binary read. The tape assign card should 

specify that the tape is an "L'' tape written at 800 bpi. 
The word count should specify one more word than the 
largest record. 

C3.0 TAPE FORMAT FOR MA-007 (UNIVERSITY OF WYOMING) 

Machine - Sigma 7« JSC does not have a Univac 1108 to Sigma 7 

conversion package. However, the University of Wyoming 
does have a CDC 3C~0 to Sigma 7 conversion package. So 
a CDC 3000 compatiole tape is supplied to the university. 

Density - 800 bpi 

Parity - odd (binary) 

Number tracks - 7 

Word length •* ^8 bits 

preceding page blank not filmed 


C2.0 TAPE FORMAT FOR MA-007 (LRC) 

- CDC 6600 

- 800 bpi 

- odd (binary) 

- 7 

- 60 bits 

~ floating point 


153 


Word type 
Number of words 

Read 

Machine 

Density 

Parity 

Number tracks 
Word length 
Word type 
Number of words 

Read 

Machine 

Density 

Parity 

Number tracks 
Word length 
Word type 
Number of words 

Read 


floating point 

2i+3 words per record, of which the last two words are 
zeros. See appendix B5.0 for the definition of the first 
24l words . 

- Use "BUFFER IN." 

C4.0 TAPE FORMAT FOR MA-048 

- CDC 3800 

- 556 bpi 

- odd (binary) 

- 7 

- 48 bits 

- floating point 

- 243 words per record, of which the last two words are 
zeros. See appendix B5.0 for the definition of the first 
24l words . 

-Use "BUFFER IN." 

C5.0 TAPE FORMAT FOR MA-059 (JSC) 


- Uni vac 1108 

- 800 bpi 

- odd (binary) 

- 7 

- 36 bits 

- floating point 

- 24l data words per record. See appendix B5-0 for the 
definition of the data words . 

- Use Fortran binary read. 
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C6.0 TAPE FORMAT FOR MA-059 (UNIVERSITY OF MICHIGAN) 


Machine 

- IBM 371/68. An IBM 360 tape is compatible with the IBM 
371/68 machine. So an IBM 360 tape is supplied to the 
university. 

Density- 

- 800 bpi 

Parity 

- odd (binary) 

Number tracks 

- 9 

Word length 

- 32 bits 

Word type 

- floating point 

Number of words 

- 243 words per record, of which the last two words are 
zeros. See appendix B5.0 for the definition of the 
first 241 words . 


C7.0 FORMAT FOR MA-083 

Machine 

- Interdata 70* An IBM 360 tape is compatible with the 
Interdata 70 machine. So an IBM 360 tape is supplied 
to the principal investigator. 

Density 

- 556 bpi 

Parity 

- odd (binary) 

Number tracks 

-7 

Word length 

- 32 bits 

Word type 

- floating point 

Number of words 

- 243 words per record, of which the last two words are 
zeros. See appendix B5-0 for the definition of the 
first 24l words . 


C8.0 FORMAT FOR MA-088 

Machine 

- Interdata 70. An IBM 360 tape is compatible with the 
Interdata 70 machine. So an IBM 360 tape is supplied 
to the principal investigator. 

Density 

- 556 bpi 

Parity 

- odd (binary) 
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Number tracks 
Word length 
Word type 
Number of words 


Machine 

Density 

Parity 

Number tracks 
Word length 
Word type 
Number of words 

Read 


Machine 

Density 

Parity 

Number tracks 
Word length 
Word type 
Number of words 


7 

32 bits 

floating point 

243 words per record, of which the last two words are 
zeros. See appendix B5.0 for the definition of the 
first 24l words . 


C9.0 TAPE FORMAT FOR MA-O89 

cdc 64oo 

800 bpi 
odd (binary) 

7 

60 bits 

floating point 

24l data words per record. See appendix B5.0 for the 
definition of the data words . 

Use Fortran binary read. The tape assign card should 
specify that the tape is an "L" tape written at 800 bpi 
The word count should specify one more word than the 
largest record. 


C10.0 TAPE FORMAT FOR MA-106 

CDC 7600. A CDC 6600 tape is compatible with the 
CDC 7600 machine. So a CDC 6600 tape is supplied to 
the principal investigator. 

800 bpi 

odd (binary) 

7 

60 bits 

floating point 

24l data words per record. See appendix B?.0 for the 
definition of the data words. 
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Read 


Machine 

Density 

Parity 

Number tracks 
Word length 
Word type 
Number of words 

Read 

C12.0 

Machine 

Density 

Parity 

Number tracks 
Word length 
Word type 
Number of words 

Read 


- Use Fortran binary read. The tape assign card should 
specify that the tape is an "L" tape written at 
800 bpi . The word count should specify one more word 
than the largest record. 


C11.0 TAPE FORMAT FOR MA-136 (DMAAC) 


- Uni vac 11 08 
-800 bpi 

- odd (binary) 

- 7 

- 3 6 bits 

- floating point 

- 2Ul data words per record. See appendix B5.0 for the 
definition of the data words . 

- Use Fortran binary read. 


TAPE FORMAT FOR MA-136 (SMITHSONIAN INSTITUTE) 


- Univac 1108 

- 800 bpi 

- odd (binary) 

-7 

- 36 bits 

- floating point 

- 2^1 data words per record. Pee annendix B5.0 for the 
definition of the data words. 

- Use Fortran binary read, 
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ASTP POSTFLIGHT EXPERIMENT SUPPORT DATA ACCURACIES 


D1.0 INTRODUCTION 


Many principal investigators (Pi's) for the Apollo-Soyuz Test Project (ASTP) 
have expressed interest in the accuracy of ephemeris, pointing, and field-of-view 
information that they will receive from the Mission Planning and Analysis Division 
(MPAD) as a part of their postflight experiment support data requirements. 

Because of the nature of the many parameters that affect the uncertainties, there 
is no straightforward method of describing these accuracies for all cases. 

Parameters that directly affect these accuracies, such as the uncertainty in the 
vehicle ephemeris, can vary radically for different mission phases. The accuracy 
of the pointing data used by an investigator may improve as his knowledge of his 
instrument rcisalinement is refined through data processing and analysis. For 
these and many other reasons, it is desirable to devise a method through which 
each investigator may determine his own expected accuracies on the basis of his 
own particular set of conditions. This report describes the error sources that con- 
tribute ephemeris, pointing, and field-of-view uncertainties , in such a manner that 
each PI may use the postflight information provided by the MPAD to compute his own 
particular expected errors. 


D2.0 SYMBOLS 


ASTP Apollo-Soyuz Test Project 

d distance 

EOM end of mission 

g acceleration of gravity at the Earth's surface 

IMU inertial measurement -unit 

MPAD Mission Planning and Analysis Division 

nav. navigation 

PI principal investigator 

rss root-sum-square 

At elapsed time 

6p position error 


- n P0..MSD 
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D3.0 ANALYSIS 


The MPAD will provide three basic classes of data to Pi’s requesting data - 
spacecraft ephemeris, spacecraft orientation, and instrument pointing. The first 
two classes are mutually exclusive (computed independently), and consequently 
their uncertainties are not related. Pointing errors, however, are a function of 
both ephemeris errors and orientation errors, as well as instrument misalinement . 

To understand this relationship, it is necessary to understand the transformations 
required to express the instrument line-of-sight definition in an Earth-fixed 
inertial system. 

The nominal instrument line-of-sight is defined by a set of alinement angles 
that relate the line-of-sight to the rigid body of the Apollo spacecraft, or the 
navigation (nav.) base coordinate system. This transformation notation matrix is 
T 

defined as [F] . The spacecraft body, in turn, can be related to the nominal 
inertial measurement unit (IMU) or reference platform coordinate system by the 

T T 

rotation defined by the matrices [d] and [G] . The G matrix relates the 
actual IMU system to the nav. base system and is computed from the spacecraft 
orientation, or gimbal angles. The D matrix relates the nominal and actual 
IMU systems and accounts for platform misalinement s and/or drift, if known. The 
nominal IMU system is related to the Earth-fixed system through a rotation matrix, 

T 

[R] . The R matrix is known as REFSMMAT. By combining these rotations, the 
nominal instrument line-of-sight can be defined in the Earth-fixed coordinate 

T 

system by utilizing the rotation matrix [FGDR] . This, then, defines the direction 
of the line-of-sight. The direction of the line-of-sight must be combined with 
the position of the instrument in Earth-fixed space to compute what the instrument 
is viewing. In practice, the spacecraft position is normally assumed as the 
instrument position, and spacecraft center of gravity/ instrument location offsets 
are ignored. The effect of this assumption is an error of 3 meters or less for 
viewing distances of at least 225 kilometers , or an angular field-of-view error of 
less than 0.001 degree. Once the relationships discussed above are understood, it 
becomes reasonably simple to determine the prime sources of error in instrument 
line-of-sight and/or field-of-view computations. 

The IMU is alined by the astronauts by observing known stars in their sextant. 
The 3a uncertainty in a nominal alinement is 0.03 degree. The platform will drift 
from thr s initial alinement. The expected 3a drift is 0.09 deg/hr. By computing 
the drift rates following each alinement and compensating for them, this uncertainty 
can be reduced to 0.02 deg/hr. All known drifts will be removed in the processing 
of instrument orientation data for ASTP experiment support. Thus, the expected 
3a uncertainty in the REFSMMAT, or R rotation, is 0.03 degree plus the uncertainty 
due to platform drift (0.02 deg/hr). The D rotation is an identity unless 
misalinement and/or IMU drift- are known and accounted for. Since the experiment 
support data will be corrected for drift, the uncertainty in the R rotation due 
to drift will be 0.02 deg/hr. For each case, the actual uncertainty due to 
drift is determined by the time interval between the latest alinement and the 
observation being performed. Thus, if observations are performed 90 minutes 
following a platform alinement, the expected IMU orientation uncertainty would be 
0.03 degree due to misalinement and 0.03 degree due to drift, or a total of Q.Oki 
degree (root-sum-square (rss)). 
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The relationship of the spacecraft or nay. "base coordinate system to the 
platform system is described by the rotation matrix ft. The expected 3a 
error in boresighting the nav. base (the mounting error in the platform relative 
to the spacecraft body) is 0.05 degree. The expected error due to inaccuracies in 
the gimbal angles is 0.01 degree per axis, or a total rss error of 0.017 degree. 
Thus, the total expected error in the G transformation is 0.053 degree. 

The transformation from the spacecraft body or nav. base system to the 
nominal instrument system described by the P transformation matrix is the 
largest error source in the line-of-sight computations . The error sources 
affecting this rotation and the corresponding error magnitudes are listed below. 


1. Instrument alinement 0.50° 

2. Mounting hole accuracy 0.25° 

3. 1-to-zero-g structural deflections 1.00° 

H. Structural thermal variations 0.33° 

Total (rss) 1.19° 


These error sources are collectively referred to as instrument misalinement . 

The rotations discussed above determine the orientation of the instrument 
line-of-sight. To determine its position in space, and thus the point being 
viewed and/or the field-of-view, the vehicle position is also required. The 
3a uncertainties in the inertial vehicle position for the ASTP experiment support 
data will range from a minimum of 600 meters to as much as 10 kilometers. 

Nominally, these errors will be about 3 kilometers, 3a. Only in cases such as 
the MA-059 data takes, where considerable command-and-service-module thrusting 
is taking place, should the uncertainties exceed 3 kilometers. 

These position uncertainties will have varying effects on the accuracy of 
line-of-sight and field-of-view calculations, depending on the distance to 
the object being viewed. For the limiting case, star observations, the target, 
for all practical purposes, is at an infinite distance and the pointing com- 
putations are unaffected by vehicle position variations. For other cases - in 
particular, for Earth observations, where the target may be as close as 225 kilometers 
- the vehicle position errors can become significant. For the case of highest 
expected uncertainty, where the maximum expected error of 10 kilometers is 
combined with a minimum target distance of 225 kilometers, the angular error in 
the co'W'nuted field-of-view could reach 2.55 degrees. 

Obviously, the effects of spacecraft position errors cannot be simply 
generalized. Each case must be examined separately, with consideration given 
to both the magnitude of the expected position error and the distance of the 
target. The WAD will provide Pi's with trajectory uncertainties for each 
segment of the mission. 
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On the basis of the discussions in the preceding paragraphs, the error 
sources that affect instrument line-of-sight • and f ield-of-view calculations can be 
summarized as follows. 


Source 

1. Platform misalinement 

2. Platform drift 

3. Nay. base misalinement 

4. Gimbal angle errors 

5. Instrument misalinement 

6. Field-of-view uncertainty 

7. Vehicle position error 

a At is the time, in hours, from 
time of interest. 

b If not corrected. 


Magnitude (degrees) 

0.03° 

( 0 . 02° /hr ) ( a A t ) 

0.05° 

0 . 02 ° 

b 1.19° 

Experiment dependent 
Experiment dependent 
latest platform alinement to the 


Sources 1, 3, and 4 are fixed and will be the same for all cases. Source 2 
varies with At and must be determined separately for each case. Source 5 is 
fixed unless the PI is able to recover a misalinement estimate from his data. 

In its processing, the MPAD will include a rotation to account for instrument 
misalinement . This rotation will nominally be an identity unless the PI is 
able to specify the misalinements prior to the processing. These will be 
required by end of mission (E0M) + 3 months at the latest, and notification 
that a nonzero misalinement is expe.cted will be required by E0M + 6 weeks. 

In the event that the misalinement error magnitude is to be corrected in the 
processing, the uncertainty estimate to replace the uncorrected uncertainty of 
1.19 degrees will have to be determined by the PI. In no cases should this 
uncertainty be reduced significantly below the 0.33 degree due to thermal 
variations. Source 6 is the uncertainty in the geometric definition of the 
instrument field-of-view. It does not affect line-of-sight calculations. Again, 
determination of the magnitude of this uncertainty is the responsibility of the PI. 
Source 7, is previously discussed, must be evaluated separately for each case and 
may vary from zero degrees, for star observations, to several degrees for local- 
vertical Earth observations. The formula for approximating the uncertainty is 


3a = sin" 1 6p/d 

where 6p is the spacecraft position uncertainty (to be provided by the MPAD) 
and d is the distance to the target. 
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The ASTP PI, "by evaluating the magnitude of the uncertainty contributed hv 
each of the seven sources noted previously and determining the root-sum-square 
of the results, can evaluate his expected pointing error for any segment of the 
mission. 
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